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Foreword ofthe Editor 


Polyethylene terephthalate (PET) is one ofthe most widely used plastics, with 
an annual primary global production of 33 million tons in 2015. 
However, only a small fraction of the used and mostly contaminated PET 
material is recycled to date. Even more, today the process of recycling leads to 
a steady deterioration in quality of the material. Excluding a costly color 
presorting process, the only known recycling process of contaminated PET is 
the chemical recycling. However, this is a costly and energy-intensive process. 
On the other hand, alkaline hydrolysis of PET would be a promising recycling 
method. In alkaline hydrolysis, PET is hydrolyzed into sodium hydroxide 
NaOH, producing the disodium salt Na2[ PTA] and monoethyleneglycol. New 
PET can be produced from the two resulting substances the para-terephtalic 
acid (PTA) and the monoethylene glycol (MEG). In this frame, the microwave- 
assisted alkaline hydrolysis of PET would add significantly to an energy 
efficient recycling process. By using dielectric heating, a significant reduction 
of the reaction time from about 3 hours to less than 10 minutes is possible 
while, at the same time, the reaction takes place at lower temperatures. This 
work presented here bases on a simple microwave applicator at laboratory 
scale that already existed. Developed by the company Gr3n, two household 
microwave ovens were modified and connected to each other. 


The aim of this work is to develop an energy-efficient modular-type and 
scalable microwave applicator to support the depolymerisation reaction, with 
a minimum capacity of 12.5 kg/h PET and a continuous operation by means of 
an Archimedean screw. To enable an efficient microwave process, all materials 
involved in the reaction are dielectrically characterized. The dielectric 
properties are systematically used in the design of the microwave applicator. 
A systematic investigation of the chemical reaction is done to make a statement 
regarding the reaction enthalpy. In order to regulate the microwave applicator, 
the modelling of the reaction kinetics is done. With the help of that a highly 
efficient industrial microwave reactor could is proposed in this work. 


Foreword of the Editor 


In frame of his work, Vasileios Ramopoulos presents several different 
innovations. In the area of the dielectric measurements an innovative, in-situ 
dual mode resonator approach is presented which, in combination with 
numerical simulation, enables a significant expansion of the measurements 
ranges for materials with different dielectric constants and loss factors. By 
using a suitable sample container, measurements above the boiling point are 
realized. In the area of calorimetric measurements, a microwave calorimeter is 
presented that enable the determination of the reaction enthalpy and the heat 
capacity. The modelling of the reaction kinetics is also presented based on the 
measured dielectric losses. Finally, a novel industrial scale microwave 
applicator at 2.45 GHz with homogeneous field distribution is proposed. 
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Kurzfassung 


Alleine im Jahr 2017 wurden 4,2 Millionen Tonnen Kunststoffabfälle von 
Europa, den USA und Kanada in Entwicklungsländer, wie zu Beispiel 
Malaysia, Indonesien und Indien exportiert, davon waren ca. 330 Kilotonnen 
PET (Polyethylenterephthalat). Meistens sind diese Zielländer nicht in der 
Lage diese Kunststoffabfälle sinnvoll weiterzuverarbeiten. Diese werden oft 
deponiert, landen im Meer oder werden durch Verbrennung energetisch 
verwertet. Dies führt zu einer nachhaltigen Beeinträchtigung des Lebensraums 
der Tiere und des Menschen. Die chemische Rückgewinnung der bei der 
Herstellung von PET verwendeten Monomere stellt die einzige Möglichkeit 
dar, PET vollständig zu recyceln. Insbesondere die alkalische Hydrolyse ist 
das einzige Recyclingverfahren, das auch für stark kontaminiertes PET wie 
Verpackungsfolien, PET-Flaschen oder Landwirtschaftsfolien geeignet ist [1]. 
Die Reaktionsprodukte der alkalischen Hydrolyse könne direkt zur Herstellung 
von neuem, reinem PET verwendet werden. Da dieses Verfahren einen 
langwierigen und energieintensiven Prozess darstellt, wird es nur selten 
eingesetzt. Im Unterschied dazu ist das Mikrowellenunterstützte Recycling 
(Depolymerisation) von PET eine vielversprechende Technologie, die die 
Reaktionszeit von ca. drei Stunden auf weniger als zehn Minuten reduziert. 
Dieses Verfahren macht das chemische PET- Recycling zum ersten Mal 
wirtschaftlich rentabel und attraktiv für die industrielle Anwendung. 


In dieser Arbeit wurde ein absolut neuer und fortschrittlicher Mikrowellen- 
Applikator zur PET- Depolymerisation mittels Mikrowellen entwickelt. Der 
Applikator wurde für den industriellen Gebrauch im ISM- Frequenzband von 
2,45 GHz ausgelegt. Bei der Entwicklung wurde Wert auf den modularen und 
skalierbaren Aufbau sowie die Energie-Effizienz des Systems gelegt. Die 
homogene Erwärmung des Reaktionsgemisches wurde mittels eines 
innovativen Wellenleiternetzwerks zur Verteilung der Mikrowellenleistung 
ermöglicht. Dieses Verteilnetzwerkkonzept ist konfigurierbar und lässt sich 
mit moderatem Aufwand auf andere Durchmesser als auch Längen des zu 
beheizenden Reaktors anpassen. Eine an das elektromagnetische Feld 
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Kurzfassung 


angepasste archimedische Schraube erlaubt einen Durchlaufbetrieb mit 
flüssigen Proben bzw. Dispersionen. Der konzipierte Applikator erlaubt einen 
Durchsatz von mehr als 100 Tonnen PET bzw. 750 Tonnen Reaktionsgemisch 
pro Jahr im Dauerbetrieb und erreicht einen Wirkungsgrad von mindestens 97 
% bezogen auf die installierte Mikrowellenleistung, wodurch er besonders 
energieeffizient ist. Je nach Anforderungen kann der Applikator modular als 
Kaskade realisiert werden, um eine längere Heizstrecke bzw. den Einsatz von 
mehreren Leistungsquellen zu ermöglichen. Um den oben genannten 
Durchsatz zu erreichen, werden beispielsweise fünf Module benötigt. 


Um den Applikator auf die Depolymerisierungsreaktion abzustimmen, ist die 
genaue Kenntnis der dielektrischen Parameter des Reaktionsgemisches 
notwendig. Dafür wurde ein innovatives Messsystem, das auf einem 
Dual-Mode Resonator basiert, entwickelt. Dieses System beruht auf der 
Störkörpermethode, die eigentlich nur für kleine Störungen, d.h. für kleine 
Proben mit geringen dielektrischen Verlusten, geeignet ist. Aufgrund des 
Dual-Mode-Ansatzes und einer vollständigen numerischen 3D- Simulation des 
inversen Problems wird die Störkörpermethode erstmals für Proben mit sehr 
hohen Verlusten eingesetzt und ermöglicht die Vermessung von Flüssigkeiten 
mit hohen Verlusten unter Druck bis zu 20 bar. Es macht den Weg frei für in- 
situ temperatur- und zeitabhängige dielektrische und kalorimetrische 
Messungen im  Mikrowellenbereich. Es lassen sich mittels 
Mikrowellenerwärmung verschiedene definierte Heizprofile in-situ 
untersuchen. Durch die kalorimetrischen Messungen ist eine Aussage 
bezüglich des exothermen oder endothermen Verhaltens der Reaktionskinetik 
möglich. Dies ist wichtig für die spätere Auslegung der installierten 
mikrowellenbasierten Leistungsquellen sowie des benötigen Kühlsystems für 
den finalen Mikrowellen-Applikator. 


Des Weiteren wurde in dieser Arbeit ein Modell zur Reaktionskinetik unter 
Verwendung des dielektrischen Verlustfaktors entwickelt. Dieses Modell 
ermöglicht die spätere Optimierung der Prozessparameter. 


Abstract 


In 2017, 4.2 million tons of plastic waste, of which approximately 330 kilotons 
were PET (polyethylene terephthalate), were exported from Europe, the United 
States and Canada to developing countries such as Malaysia, Indonesia and 
India. Most of these target countries of plastic waste are unable to process them 
sensibly. These plastic wastes are often landfilled, end up in the ocean or are 
burned for energy. This leads to a continuous disturbance of the habitat of 
animals and humans. Chemical recycling of PET is the only way to fully 
recover the monomers used in the manufacture of PET. In particular, the 
alkaline hydrolysis is the only recycling process that is also suitable for highly 
contaminated PET such as packaging films, PET bottles or agricultural films 
[1]. The reaction products of the alkaline hydrolysis can be used directly for 
the production of new, virgin PET. Since this procedure is a lengthy and 
energy-intensive process, it is rarely used. In contrast to this, microwave- 
assisted recycling (depolymerization) of PET is a promising technology that 
reduces the processing time from around three hours to less than ten minutes. 
This process makes chemical PET recycling economically profitable and 
attractive for the industrial application. 


In this work, a new microwave applicator for PET depolymerization on an 
industrial scale is demonstrated. The applicator is designed for industrial use 
on the ISM band frequency of 2.45 GHz. During the development, the focus is 
put on the modular and scalable design as well as the energy efficiency of the 
system. The homogeneous heating of the reaction mixture is made possible by 
an innovative microwave power distribution network. This distribution 
network concept is configurable. Both the diameter and the length can be 
adapted with little effort, depending on the dimensions of the reactor to be 
heated. An Archimedean screw adapted to the microwave field allows 
continuous operation with liquid samples and dispersions. The applicator 
enables a throughput of more than 100 tons of PET or 750 tons of reaction 
mixture per year in continuous operation. The power absorbed in the sample is 
more than 97 % compared to the installed microwave power, which makes it 


Abstract 


very energy-efficient. Depending on the requirements, the applicator can be 
implemented as a modular cascade to enable a longer heating section or the use 
of several power sources. For example, five modules are required to achieve 
the throughput mentioned above. 


To adapt the applicator to the depolymerization reaction, a precise knowledge 
ofthe dielectric parameters of the reaction mixture is necessary. A dual-mode 
cylindrical cavity is developed for this purpose. This system relies on the cavity 
perturbation method, which is actually only suitable for small perturbations, 
i.e. small samples with low dielectric losses. The dual-mode approach in 
combination with a full 3D numerical solution of the inverse problem allows 
the measurements of liquids with high losses under pressures up to 20 bar. At 
the same time, it allows in-situ temperature and time dependent dielectric and 
calorimetric measurements at 2.45 GHz. By means of microwave heating 
different defined heating profiles can be investigated in-situ. The calorimetric 
measurements provide information about the exothermic or endothermic 
behavior of the reaction kinetics. This is important for the later design of the 
installed microwave power sources as well as for the required cooling of the 
microwave final applicator. 


Furthermore, in this thesis a model for the reaction kinetics using the dielectric 
loss tangent was developed. This model allows the optimization of the process 
parameters. 
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1 Introduction 


1.1 Dielectric heating using microwave power 


Today, the use of microwaves for dielectric heating finds increasing 
application in heating of materials such as food [2], in microwave assisted 
chemistry [3], [4], [5], [6] and in drying of materials [7]. The advantage of 
microwave heating originates from the fact that the heat is generated inside the 
material volume, so the heating is quick, direct and energy efficient. In 
particular, materials characterized by low thermal conductivity such as 
plastics, glass, ceramics and powders can be heated much faster than with 
conventional methods. Furthermore, microwaves enable the selective heating 
of mixtures of materials, having different dielectric loss factor. This allows 
process specific temperature profiles that are impossible to realize in 
conventional thermal furnaces. Another advantage of microwaves is that 
workflows at very high temperatures can be realized in cold ovens. 


In order to develop energy efficient microwave systems and processes, the 
detailed knowledge of the dielectric properties of the processed materials at the 
center heating frequency is a prerequisite. The lack of data on temperature 
dependent dielectric properties possess barriers for further development as they 
can be rarely found in literature for the temperature and frequency of interest. 
Even if such data can be found, very often the reported material differs from 
the specific material of interest. Therefore, the accurate simulation of processes 
and microwave applicator demands the precise dielectric characterization of 
all involved materials. In case the materials exposed to microwave heating pass 
through irreversible changes, for example in chemical reactions with 
consecutive change in dielectric properties, the knowledge about the time and 
temperature history of material’s permittivity is essential. 


Microwave assisted chemistry is reported to provide numerous advantages, 
such as remarkable acceleration of the reaction, higher yields under milder 
reaction conditions and higher product purities [8]. These advantages may be 
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explained as a combination of thermal effects, such as higher heating rates, 
superheating (hot spots) and selective absorption ofthe microwave power from 
different substances. 


1.2 Motivation 


Polyethylene terephthalate, PET, is one of the engineering plastics which is 
widely used as it is showing excellent thermal and mechanical properties, 
chemical resistance, and transparency. By far the largest application is in the 
textile industry. Enormous quantities of this material are also consumed in the 
manufacture of PET films, used in 3D printing, injection molding, food 
packaging, and especially for the production of water and soft drink bottles. 
PET is the sixth-most used plastic with a total primary global production of 33 
million tons in 2015 [9]. 


The usual durability of polymers such us PET is, on the one side, desirable for 
its application. On the other side, it is a threat to the environment when these 
are disposed after use. Different to the natural polymers, the synthetic polymers 
are not easily degraded under normal outdoor conditions. They are collected in 
landfills as waste products [10]. At the same time, discarded polymers 
represent a huge waste due to the amounts of energy and resources embodied 
in them [10]. Therefore, the recycling of PET becomes inevitable. The 
recycling of PET can be realized in many ways: (a) primary recycling refers to 
the mechanical reprocessing of uncontaminated polymer products into 
products with equivalent properties; (b) secondary mechanical recycling where 
the polymer is separated from the various contaminants and can be reprocessed 
into products with inferior properties; (c) chemical recycling refers to the full 
recovery of the monomers and (d) energy recovery [11], [12]. The chemical 
recycling of PET is the only ‘true’ recycling technique, which does not lead to 
a degradation and it is suitable for contaminated polymer products [12]. 


The chemical recycling is based on processes such as hydrolysis, methanolysis, 
glycolysis, ammonolysis and aminolysis [13]. Depending on the reaction 
chosen, different reaction products are obtained as shown in Fig. 1.1. Until 
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today, only the glycolysis and methanolysis can be commercially applied for 
the depolymerization of PET [10]. 


Chemical 
recycling 
of PET 


Other 
Methanolysis: Hydrolysis: Processes: 
DMT + MEG PTA+ MEG Aminolysis 

Ammonolysis 


Glycolysis: 
BHET + 
Oligomer 


Neutral Alkaline 
hydrolysis hydrolysis 


Acid hydrolysis 


Figure 1.1: Possible chemical reactions for the depolymerization of PET and the 
corresponding reaction products 


The hydrolysis is the only recycling method, which results in base products 
such as para-terephthalic acid (PTA) and mono-ethylene glycol (MEG), which 
can be used to produce virgin PET again. 


Depending on the catalysts used, a distinction is made between acidic, alkaline 
and neutral hydrolysis. The acid hydrolysis is mostly realized in the presence 
of sulphuric acid, nitric or phosphoric acid. This process is very costly due to 
the demand of the further recycling of large H,SO, amounts and the separation 
of MEG from sulphuric acid [10]. The neutral hydrolysis is realized with water 
or steam as catalyst. The main drawback from this reaction is that the products 
show a low purity [10]. The alkaline hydrolysis is performed with the use of 
alkaline solution of NaOH or KOH [10]. The alkaline hydrolysis may be 
carried out with highly contaminated PET and it is much cheaper and faster 
compared to methanolysis. 
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Based on above stated, the alkaline hydrolysis of PET is the only interesting 
recycling method of contaminated post-consumer PET. However, the long 
reaction time in combination with the high-energy consumption is the main 
concern of this recycling technique. Because of that, PET material used today 
is only reused to approximately 10 % but not really recycled or depolymerized. 
This motivates the investigation of alternative procedures to make the alkaline 
hydrolysis of PET more time- and energy-efficient. It leads to the microwave- 
assisted depolymerization of PET that provides a significant acceleration of 
the reaction process time from about 180 minutes by means of conventional 
heating, down to less than 10 minutes using microwave heating. Additionally, 
the microwave assisted depolymerization takes place at lower temperatures 
and therefore under lower pressure conditions because of the volumetric and 
selective heating characteristics of microwaves [12]. These facts enable an 
easier upscaling of a future reactor. 


To enable a microwave-assisted reaction process, a suitable process chamber 
is required. In general, those microwave applicators are still based on the 
design of Raytheon developed in 1965 [14]. That microwave oven consists of 
a magnetron [15] as power source and a rectangular cavity of some 40 cm 
height, 50 cm width and 30 cm depth [14]. The main drawback for industrial 
application is the difficulty to enable a well-defined and uniform heating of 
materials, because of the existing cold and hot spots in any microwave oven. 
This problem of non-uniform heating until now is one of the main obstacles 
for the commercialization of microwave heating in industrial applications [16]. 
The main factors, which contribute to this lack of uniform heating, are the 
geometry of the microwave resonator, the finite penetration depth of the 
microwaves, the geometry of the material sample, the dielectric properties of 
the processed material and the feed point of the microwaves compared to the 
position of the processed material [17]. These factors make it also impossible 
to get a universal microwave oven for the processing of any dielectric material 
and geometry. 


This motivated the development of various lab scale and industrial scale 
application specific microwave ovens with the main target to offer more 
homogeneous field patterns [16], [18], [19] and [20]. In [20] a novel industrial 
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microwave system is presented, which is based on a patented hexagonal cavity. 
This geometry provides a significantly increased homogeneity in a large 
applicator volume [21]. In [22], a ‘Coaxial Wrap-Around TE), Mode 
Converter’ cavity is based on a TE), waveguide in a circular orientation. The 
power is coupled into the actual heating cavity through nine radial slots (see 
Fig. 3 in [22]). Further industrial microwave applicators for both batch and 
continuous operation are shown in [23], [24] and [25]. In [18] various 
microwaves ovens are presented, which support the microwave-assisted 
organic synthesis. Those ovens are based on single and multimode cavities and 
are partially equipped with mode stirrers, special temperature sensors and 
pressure vessels to enable well-defined and reproducible heating conditions. 
These systems are mainly targeted to lab-scale operation. 


In the area of microwave assisted PET recycling based on single mode cavities 
at lab-scale several investigations are on-going [26], [27] and [28]. Most of 
them are based on commercial single mode microwave reactors developed 
from CEM corporation [29]. Those investigations were carried out under 
atmospheric pressure and under higher pressure with samples between 10 ml 
and 100 ml. However, the main problem for the design of an industrial scale 
applicator is the restricted penetration depth (in the range of some millimeters) 
in the reaction mixture and the demand for a continuous operation to meet the 
economic needs. 


1.3 Objectives and structure of this work 


The focus of the current research is to develop an efficient industrial scale 
microwave reactor with a PET material throughput of more than 157 t/a and 
accordingly a total material throughput from about 1200 t/a to support the 
microwave assisted alkaline hydrolysis of PET. In the patents [30] and [31] a 
new method, the corresponding apparatus and scenario for a continuous 
heating of a high-loss reaction mixture with an electromagnetic field are 
presented. The design relies on a metallic Archimedean screw placed in a 
microwave transparent tube, which transports the full mixture through the 
microwave cavity. The present research is relying on the patents [30], [31] and 
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on a new microwave applicator design, which meets the needs of an industrial 
scale process of depolymerization of PET. 


Detailed knowledge of the dielectric properties of the material parameters is 
essential for successful system and process design. Therefore, a system for 
in-situ monitoring of the chemical reaction under microwave heating 
conditions must be developed. This system shall allow the reconstruction of 
the temperature dependent dielectric properties of the reaction mixture 
consisting of mono-ethylene glycol, sodium hydroxide, water and PET in the 
ISM band at 2.45 GHz. In addition, this system must allow in-situ temperature 
and time dependent dielectric measurements under microwave heating 
conditions. It must be taken into account that the reaction takes place at 
temperatures above the boiling point and the system must accordingly be able 
to keep the reaction mixture under pressure. In addition to the dielectric 
measurements, a calorimetric measurement setup under microwave heating 
conditions must be developed. These measurements will be used to identify 
the exothermic or endothermic characteristics of the proposed reaction. 


Secondly, based on the measured dielectric properties, a modular and scalable 
applicator on industrial scale has to be developed. This applicator must allow 
continuous operation using a metallic Archimedean screw within a glass tube, 
transparent to microwaves. In this context, the screw geometry has to be 
optimized with respect to the desired production rate and the overall efficiency 
of the microwave power. Furthermore, the proposed microwave cavity must 
ensure a homogeneous field distribution in the axial and azimuthal direction of 
the reactor. 


In addition, a process optimization scenario must be introduced. Thereby the 
measured temperature and time-dependent dielectric properties have to be used 
for development of an appropriate reaction kinetics model. In this way, the 
different reaction scenarios can be evaluated to obtain an optimal reaction 
process. 


Based on the listed objectives, chapter 2 describes the theory of dielectric 
properties, applied dielectric mixing rules and the alkaline hydrolysis of PET. 
The developed dielectric and calorimetric measurements setup and the 
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corresponding calibration procedure are presented in chapter 3. The measured 
dielectric properties of the individual reactants and of the full reaction mixture 
are described in chapter 4. Chapter 5 provides an overview of the developed 
model of the reaction kinetics based on the measured dielectric properties. The 
design of the developed industrial scale microwave cavity and experimental 
results with the developed reactor are shown in chapter 6. Finally, a conclusion 
and outlook are given in chapter 7. 


2 Fundamentals 


In this chapter, the microwave-assisted depolymerization of PET based on high 
temperature alkaline hydrolysis are introduced in detail. First, the principles of 
microwave heating are addressed. Then the different measurements methods 
for the dielectric characterization of materials and the alkaline hydrolysis of 
PET are explained. 


2.1 Definition of dielectric properties 


In general, electromagnetic heating refers to three different heating 
mechanisms, which are described by the material-specific parameters. The 
interaction of materials with electromagnetic waves is characterized by their 
dielectric (permittivity), magnetic (permeability) and conductive (electric 
conductivity) properties. Usually these material parameters depend on the 
frequency of the electromagnetic field, the temperature of the material and in 
case of anisotropic properties, as well on the orientation of the material versus 
the direction of the electric field [32]. Depending on the dominant material 
parameters, the electromagnetic heating may be realized through dielectric 
heating and/or magnetic heating and/or conductive heating. 


Dielectric heating relies on the ability of the material to absorb energy when 
an external alternating electric field is applied. Such a material is called 
‘dielectric material’ [32]. A dielectric material is characterized by its complex 
dielectric constant €. The complex dielectric constant € consists of a real part 
e’, which represents the energy storage term, and an imaginary part e” which 
represents the dielectric loss term. 


E€ = €' — je" = €o (€p — jer) (2.1) 
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where e, = 8.8542 - 1071? C?/(N m?) is the vacuum permittivity and €; 
and €, are the relative dielectric constant and relative dielectric loss factor, 
respectively. 


Magnetic heating relies on the ability of a magnetic material to absorb energy 
when an external alternating magnetic field is applied. A magnetic material is 
characterized by its permeability u. The complex permeability u as well consist 
of a real part u’ that represents the energy storage term and an imaginary part 
u” that represents the magnetic energy loss term [33]. 


w= — ju" = Wo (Uy — jur) (2.2) 


where uo = 47-1077 H/m is the vacuum permeability and w. and uf are the 
relative magnetic constant and relative magnetic loss factor, respectively. 


Conductive heating relies on the current flow, which takes place in a 
conductive material when an external electric field is applied. A conductive 
material is characterized by its electronic and/or ionic conductivity ø. Since 
many materials or mixtures of materials comprise dielectric as well as 
conductive losses, often the effective conductivity o,¢¢ is given as the sum of 
dielectric and conductive losses by the following formula [32]: 


Ceff =T+E"W (2.3) 
where w is the applied field angular frequency. 


A dielectric material may experience a variety of polarization mechanisms and 
loss mechanisms. These parameters strongly depend on the material 
composition as well as the operating temperature and frequency range [33]. 
The loss and polarization mechanisms may be classified into ionic conduction, 
dipolar relaxation, atomic polarization and electronic polarization (see Fig. 
2.1). The sum of all mechanisms results in the overall permittivity. In the 
microwave frequency range, the ionic conductivity and the dipolar relaxation 
mainly determine the loss behavior of a dielectric material. The effective 
material losses can be described as the sum of dielectric losses €” and 
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conductive losses ø. Depending on the materials composition the conductivity 
is based on electronic or ionic conductivity. 


| tv 


> Dipolar 
€, Rotational) 


ig 


Atomic Electronic 


Figure 2.1: Dielectric permittivity over a wide frequency range [32] 


The ratio of €; to €, is called dielectric loss tangent tan 6 [32]. 


A 
r 


tan ô = (2.4) 


er 


The absorbed power as a function ofthe applied electric field is determined by 
the following formula [17]. 


Days | we' tan ô E? dV (2.5) 
Vs 


V, is the volume of the dielectric material and E is the amplitude of the 
electrical field. If the electromagnetic wave propagates in a lossy dielectric 
material, its amplitude progressively attenuates in the propagation direction. 
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The distance at which the power of incident wave attenuates by l/e, 
approximately 37 %, is called penetration depth D,. Assuming that u = 0, it 
can be approximated as follows [34]: 


a 
pee Z [VIF tan®5 - 1] 2 (2.6) 


P Anf Jet 


where Cy = 1// €o" Ho = 3: 108 m/s is the speed of light in vacuum. The 
higher the losses of the material the lower is the penetration depth. 
Accordingly, the dielectric properties of the material (e, and e€;’) have 
significant influence on the field distribution in the material. 


2.1.1 Dielectric mixing rules 


In general, the effective dielectric properties of a material mixture can be 
estimated from the dielectric properties of each constituent. For this purpose, 
several dielectric mixing rules can be applied [35]. 


The most common dielectric mixing rules are the Maxwell-Garnett (MG) 
formula and the Lichteneker-Rothers (LR) formula [35], [36]. In the case of 
the MG rule, the mixture to be analyzed consists ofthe background material in 
which spherical inclusions are embedded. The MG formula is given as: 


Ereff — Ere Eri — Ere 
Zane SHEETS (2.7) 
Ereff + 2Er e Eri + 2Er e 
where Ereffs Ere and €p; are the effective permittivity of the mixture, 
background material and inclusion, respectively. @; is the fractional 
volume of the inclusion. However, it must be mentioned, that the equations are 
based on the assumption that all inclusions are spherical [35]. The 
LR formula can be applied in a wider range than the MG formula to both 
isotropic and anisotropic mixtures and it is expressed as: 
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' 1/s 
Ener = (di; + A - piee) (2.8) 


where ¢ is the LR exponential factor and is specified by the composites micro- 
geometry [36]. 


2.2 Dielectric measurement methods at 
microwave frequencies 


During recent years, a variety of methods for measuring the dielectric 
parameters of materials have been developed [37], [38]. Usually, the choice of 
the most appropriate measurements technique depends on factors such as the 
range of the expected dielectric properties, the operating frequency, the 
required measurement accuracy, the sample geometry, and finally, the 
aggregate state of the material under test. The methods can be grouped in the 
reflection method, the transmission-reflection method and the resonance 
method as shown in Table 2.1 [33]. 


2.2.1 Reflection method 


This method is based on the measurement of the reflection of an 
electromagnetic wave entering a dielectric material A very common 
implementation of this method is the so-called coaxial probe measurement. 
The key component of this measurement method is an open-ended coaxial 
waveguide (Fig. 2.2) [32]. The material is measured either by inserting the 
probe into a liquid sample or by touching the flat surface of a liquid or by half- 
solid material. A perfect contact between the sample and the probe without any 
air gap is necessary and it is necessary, that the area and the thickness of the 
sample are sufficiently large, so that the complete fringe field of the probe is 
contained. This method enables the measurement of liquid and semi-liquid 
materials over a very broad frequency range with a single probe. It is applicable 
for medium- and high-loss dielectric materials (tan ô > 0.05). 
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2.2.2 Transmission-Reflection method 


Another method is based on the measurement of the complex transmission 
coefficient (S,,) and the complex reflection coefficient (S41) from which the 
complex material permittivity can be calculated. A sample partially fills a 
transmission line, which is usually a section of a rectangular or coaxial 
waveguide [39]. 


Table 2.1: High frequency dielectric measurement methods and their pros and cons 


Method Se Advantages Disadvantages 
Measured par. 
Coaxial Probe 
a Sample surface 
S s Broadband, quality, tempera- 
8 At Flexible sample | ture restrictions 
rm Er 
2 geometry through the 
coaxial probe 

Waveguide 
S| els $5 Complex sample 
5| == Int Broadband geometry 
Q| coaxial waveguide Er, Uy i 
z= preparation 
O 
7 
x waveguide 
2 Free space 
n 
E Sutanleorhigh | 7 SMUT 
= PR: temperatures atile 
S au reflections 
= ` 7 fr easy sample between the 

Se preparation antennas 
Antennas 
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Cavity 
High precision, |Single frequency, 
511,521 microwave heat- |high demands on 
Er, Ur ing in the cavity |sample geometry, 
3 possible small permittivity 
s 
7 
© 
m Sample Be 
115851 ; en complex sample 
<i > a High precision Br pas. 
Open resonator A <« MUT 


mitted signal, the complex permittivity of the sample can be estimated 
[32].This method enables the dielectric characterization of medium and high 
loss dielectric materials over the full frequency band of the waveguide. 
Furthermore, it allows the measurement of magnetic material properties with 
the same measurement setup. Similarly, the measurement takes place in the 
free-space method. Here antennas are used to focus microwave energy onto or 
through a sheet of material. This method is contactless and can be used on 
materials that are exposed to high temperatures [32]. 


2.2.3 Resonance method 


The resonance method is the most precise analysis for determining the complex 
permittivity and permeability of materials. However, there are strong 
limitations to the measurement technique given by the accuracy of the sample 
geometry and the maximum sample dielectric and magnetic losses. Different 
types of resonators such as cavity resonators [40], split-cylinder resonators [41] 
and Fabry-Perot resonators [42] can be used. This method allows the dielectric 
characterization of samples at the individual frequencies of the existing modes 
in the resonator. The most common resonance method is the cavity 
perturbation method, which is explained in more details as follows. 
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2.2.4 Microwave cavity 


A microwave cavity can be described as a hollow structure enclosed by 
metallic walls. In simple terms, it can be expressed as a waveguide that is short- 
circuited at both ends. The cavity may be either hollow or filled with a 
dielectric material. Standing waves are formed within the cavity at specific 
resonance frequencies. Each resonance frequency corresponds to a certain field 
distribution. This field distribution is called resonance mode. Depending on 
whether the electric field or the magnetic field is perpendicular to the third 
dimension of the cavity, a distinction is made between TE,,,; modes and 
TMnmi modes [43]. In the case of a coaxial cavity, a TEM,,,) mode can also 
be formed in which both the electric and the magnetic field are perpendicular 
to the third dimension of the cavity. The indices n, m and | correspond to the 
number of field variations along the three dimensions of the cavity [43]. 
Sometimes not all boundary conditions can be fulfilled by pure TE- or TM- 
modes. In such cases, linear combinations of TE- and TM- modes can provide 
a complete and general solution, the so-called hybrid modes [43]. 


The most common microwave cavities are based on a rectangular or cylindrical 
waveguide short-circuited at both ends. The resonance frequency of a 
rectangular cavity is determined using the following equation [43]: 


Co NTN? MTN? In\? 
fam = an + (=) + (=) (2.9) 


where Cy is the speed of light in vacuum. a, b and d correspond to the three 
dimensions of the rectangular cavity. 


Similarly, the resonance frequency of a cylindrical cavity is determined using 
the following equations [43]: 
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fami = (2.10) 


Where Pym is the root of the first kind Bessel function Jpm and Pym is the root 
of the derivative of the first kind Bessel function Jpm. a and d corresponds to 
the radius and length of the cylindrical cavity. 


Fig. 2.2 shows the mode chart of cylindrical cavities [43]. As can be seen for 
specific ratios between radius and length of the cavity more than one mode can 
be excited simultaneously at exactly the same frequency. For example, at the 
red circle marks an area where the diameter is equal to the cavity length, the 
TMo 9-mode and the TE,,,-mode can be excited simultaneously. 
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Figure 2.2: Mode chart of a cylindrical cavity [43] 
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The quality factor of a cavity with a resonance frequency wọ is defined as the 
ratio of the average energy stored in the cavity to the energy loss during one 
oscillation period [43]: 


Qo = op ete Q.11) 
loss 

where Wn and W, corresponds to the magnetic and electric energy stored and 
Pjoss Corresponds to the resonator losses. The resonator losses are the sum of 
the conducting wall losses, dielectric losses, radiation losses and coupling 
losses. Within a cavity with a very large quality factor, a small power signal 
results in a very large electric and magnetic field strength near the resonant 
frequency [44]. Such cavities are used as microwave ovens, filters etc. 
The power losses at the conducting walls of a cavity are defined as follows: 


R 
P. = =| |H?|ds (2.12) 
2 5 


where R, = ‚/wuo/20 is the surface resistivity of the metal walls and H, is the 
tangential magnetic field on the resonator walls. 


The Q-factor of a rectangular cavity with lossy walls, filled with a lossless 
dielectric material and without consideration of radiation losses and coupling 
is determined as follows [43]: 


2W We 
Qe P. 
_ kĉabdnņ 1 
— An?R, [(2l2ab bd La d (2.13) 
(Ge ) + (22) + (az) +) 
_ (kad)*bn 1 


2n?R, 2l?a3b + 2bd3 + l2a3d + ad? 


a, b and d correspond to the physical dimensions of the rectangular cavity 
resonator, 7 is the wave impedance in the medium, l is the number of half 
waves in the third dimension of the rectangular cavity and k is the wave 
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number. The Q - factor of a cylindrical cavity with lossy walls, filled with a 
lossless dielectric material and without consideration of radiation losses and 
coupling walls for a TE„mı-mode is determined as follows [43]: 


g= WoW, 
Go P, 
_ a a a 
ee en 


a and d are the cavity radius and length, respectively. Phm is the zero position 
of the Bessel function of first kind, k is the propagation constant and £ is the 
phase constant. 


Similar, the Q-factor of the TM mı-modes is given by [45]: 


211/2 
| (2.15) 
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Figure 2.3: Normalized unloaded quality factor of a cylindrical resonator for 
various modes [43] 


Fig. 2.3 shows the quality factor of an empty cylindrical cavity depending on 
the mode and the cavity geometry without consideration of coupling and 
radiation losses. As can be seen for specific ratios between radius and length 
of the cavity, the Q. of existing modes in the resonator can be influenced. 


Cavity perturbation method 


The perturbation of the field in the microwave cavity can be realized by the 
modification of the cavity volume or the introduction of a small sample of a 
dielectric material. The latter is applied for the dielectric characterization of 
materials. The material perturbation method is based on the comparative 
analysis of electromagnetic properties and is based on the assumption that the 
electromagnetic field in the cavity with and without sample remains more or 
less constant. The determination of the sample parameters is realized by the 
measurement of the frequency shift and the quality factor change, which is 
generated through the perturbation of the electromagnetic field. 
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The change of the complex angular frequency Qm of an empty cavity by 
introduction of a material sample is given by the following equation [33]. 


Am = Omo 


mo Jy, Eò By + Ho Hy) dV 


— Jy, (€r — DES E, av 


1 


Q 


(2.16) 


Am, and Am, are the complex resonance frequencies for the empty and 
perturbed cavity, respectively. e, is the complex relative permittivity Ho, H, 
and Ey, Eı are the complex magnetic fields and complex electric fields, 
respectively before and after the introduction of the sample. The denominator 
represents the integral over all the cavity volume V., while the numerator is the 
integral over the sample volume Vz. 


At the resonance frequency the stored magnetic energy W,, in the cavity is 
equal to the stored electric energy W,. Under this assumption, the Eq. 2.16 can 
be simplified as follows. 


— f, (er — DES E, dV 
u nal len an (2.17) 
Qing 2 f, (Ei E,) dV 


The complex resonance frequency Qm of a resonator is connected to the 
resonance frequency f and the Q-factor of the resonator as follows [33]: 


Am = Qm, + jm; (2.18) 
Am =20f (2.19) 

Am 
Q=; 2.20 
Zils, (2.20) 


If the perturbation of the field is small, Qm, ~ Qm, and On; >> Qm, then 


according to the cavity perturbation theory, the Eq. 2.17 can be separated in 
the real and imaginary parts. The dielectric properties of the material under test 
(MUT) are calculated as follows [33]. 
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‚_1Ih-f)% 
NL (2.22) 
nn | 


Qo and Q, represent the corresponding quality factors. V, and V, are the 
volumes of the cavity and the material sample, respectively. A and B are 
specific factors based on the geometry of the cavity and the cavity mode 
excited. For a small material perturbation (Ey = E,) A and B are constant. 


Temperature dependent dielectric measurements 


The reconstruction of the dielectric properties of several materials has been 
investigated with various methods in several publications [46], [47], [48], [49] 
and [50]. In parallel a variety of methods for measuring the dielectric 
parameters of materials over predefined heating scenarios have been developed 
during recent years [37], [38]. As it is well known, the resonant method, such 
as the cavity perturbation method is well suited for dielectric material with low 
and moderate dielectric properties [33], while the reflection method, such as 
the coaxial probe, and the transmission/ reflection method are suitable for the 
measurement of dielectric materials with medium and high losses [51]. In 
particular, the resonance method is well suited for microwave heating of lossy 
as well as low loss dielectric materials due to the high electric field amplitude 
in the cavity at resonance frequency. Of course, a sufficient source of 
microwave power must be available for that purpose. 


Measurements on dielectric materials in a broad temperature range based on 
the cavity perturbation method are presented in [38]. This experimental setup 
is based on a rectangular TE, 94 cavity. Here the sample is heated stepwise in a 
conventional oven up to 1200 °C and afterwards is moved into the cavity for 
measurements. The dielectric properties were reconstructed using 3D full wave 
simulations by solving the inverse problem. This means that the material 
parameters of MUT are varied in expected permittivity range and the Q-factor 
and resonance frequency change is simulated. 
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Couderc et al. [52] developed a dual-mode cylindrical cavity based on a 
TMgoı0o-mode and a TE,,,-mode. Here the TMgı0-mode is used as a heating 
mode at 2.45 GHz and the TE,,,-mode is applied for dielectric measurements 
at about 3 GHz. Both modes are individually coupled into the cavity with 
minimum cross-coupling. Another study based on a dual-mode coaxial cavity 
for high-temperature measurements of small samples is presented in [53]. In 
both cases, microwave tuners are required during heating to adjust the 
resonance frequency of the loaded cavity to the fixed frequency of the 
magnetron power source in order to allow sample heating. 


Nikawa et al. already developed a simplified in-situ technique for 
simultaneously heating and measuring of dielectric parameters using a 
2.45 GHz TMo;9-mode cavity [54]. For this purpose, the output power from 
the vector network analyzer (NWA) is amplified up to 100 Watt. Another 
in-situ test set is presented in [37]. This work demonstrates that it is possible 
to create a system, which enables simultaneous heating and measuring using 
the same microwave source. This system is based on a cylindrical TE,,,-mode 
cavity and on a tunable high-power microwave amplifier. Nevertheless, a 
challenge of this measurement approach is the need of strong coupling for 
sufficient heating of the material in the cavity. But a strong coupling has a 
significant influence on the measurement accuracy and so far, inhibited 
accurate measurements of the dielectric loss factor, because the coupling is 
assumed as negligible in the equations for the perturbation theory. To increase 
the efficiency and control quality, Catala et al. extended the above-described 
approach in his work [55]. A dual-mode cylindrical cavity operating at both 
TE,,, - and TMy;9-modes was developed. In the design, the TMg;9-mode was 
used for heating of the sample while dielectric measurements were performed 
using the TE,,,-mode. However, this measurement setup is restricted to 
dielectric materials with low and medium losses since the reconstruction of the 
dielectric properties is based on an analytical model. 


As already stated, the classical cavity perturbation method is valid for 
measurements of low-loss dielectrics where the perturbation of the resonance 
is sufficiently small. High-loss materials, when being placed in the cavity, 
result in a deterioration of the transmission resonance peak in such a way that 
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the classical cavity perturbation method is no longer valid. The deterioration 
of the resonance peak, however, also depends on the particular field mode 
excited. A dual-mode cavity approach going along with a calibration procedure 
based on full-wave simulation may help to extend the measurement range 
towards high loss materials and increase the accuracy of the measurements. 


2.3 Alkaline hydrolysis of PET 


In the depolymerization reaction, the PET is hydrolyzed in sodium hydroxide 
(NaOH) that yields to the disodium salt, Na2[PTA] and mono-ethylene glycol 
(MEG) [56]: 


[PET] + 2NaOH > Na, [PTA] + [MEG] (2.23) 
O O 
Il Il 
O—C— —C— OCH,CH;--.+ 2n NaOH 
PET n 


n nacoc—{_\— COONa + n C,H,O, 


The reaction mixture consists of a sodium hydroxide solution and PET 
particles. The reaction takes place at the surface of the PET flakes and therefore 
the reaction is directly influenced by the particle size. 


According to various studies, the alkaline hydrolysis of polyethylene 
terephthalate is described as both slightly endothermic [57] and slightly 
exothermic reaction [58]. According to [57], during the hydrolysis of PET an 
additional thermal energy input (endotherm) of about 5 kJ/mol is required. In 
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contrast, according to [58], an energy of about 70 kJ/mol is released (exotherm) 
during the hydrolysis of PET. This makes the further investigation of the 
reaction enthalpy particularly interesting. 


Calorimetric measurements 


In general, calorimetry is a tool to measure the heat transfer in a material 
associated with changes in the material due to a chemical reaction, physical 
change, or phase transition. In particular, by chemical reaction, calorimetry 
enables the investigation of the power required to follow the reaction 
(endothermic reaction) or the released power during the reaction (exothermic 
reaction). The most common methods of calorimetry are differential scanning 
calorimetry [59], isothermal titration calorimetry [60] and isothermal micro 
calorimetry [61]. For all of these conventional calorimeters a thermally isolated 
system is required in order to get an accurate information. 


In addition to the above mentioned calorimetry methods, there exists also 
various investigations about thermal analysis setups under microwave heating 
conditions. In this case, quantitative results are difficult to obtain because the 
thermal analysis always assumes that the sample is in thermal equilibrium with 
its environment. But this is never the case with microwave heating. Parkes et 
al. [62] have developed a thermal analysis system based on a single-mode 
rectangular cavity and a narrow band high frequency power source (see Fig. 
2.5). Since the applied power source is not frequency-tunable, the proposed 
system enables the microwave heating of the sample based on complex tuning 
elements (e.g. variable coupling iris, short plate tuner and a stab tuner). In this 
work a qualitative thermal analysis has been carried out through the sample 
temperature monitoring and the proposed applied power. Nesbitt et al [37] has 
developed a more accurate microwave calorimeter based on TE,,,-mode 
cylindrical cavity and a tunable solid-state amplifier. This system enables 
accurate power measurements with the help of a power meter. Here the 
incident and the reflected microwave power had been continuously monitored. 
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300 watt power sliding 
generator short 


thermocouple 


Figure 2.4: Calorimetric measurements under microwave heating conditions, proposed by 
Parker et al. [62] 
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Microwave chemistry shows increasing applications in the last years. It enables 
a faster and more energy efficient reaction [63]. For the successful design of 
microwave applicators, the a-priori knowledge of the temperature dependent 
dielectric properties of the involved materials is required. Furthermore, the 
detailed knowledge of the reaction kinetic and power requirements is 
important. In this chapter, a novel dielectric and calorimetric measurement 
system is presented. 


3.1 In-situ, high pressure and high precision 
measurement setup 


3.1.1 Design methodology 


The equations 2.15 and 2.16 are derived without taking into account the 
coupling ports, sample holder and perturbation of the electromagnetic field in 
the cavity when introducing the material under test. However, for strong 
material perturbation (Ey # E,) the changes of the electric field and coupling 
losses have to be taken into account. Those effects may be accounted in the 
factors A and B in equations 2.15 and 2.16, which become dependent in that 
case on the dielectric properties of the sample. Moreover, for strong 
perturbations the separation of the complex frequency into real and imaginary 
parts is no longer valid. A calibration procedure is necessary based on full wave 
numerical simulations [38]. The frequency shift Af of the resonance and/or the 
change in quality factor saturate with increasing dielectric constant and/or loss 
factor. However, this is dependent on the wave mode excited in the cavity. In 
Fig. 3.1 the electric fields for the first two modes (TE,,; and TMo,9) in a 
cylindrical cavity can be seen. Both modes are often used for the dielectric 
characterization of materials, because they have an electric field maximum in 
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the center of the cavity. Fig. 3.1 also shows the electric field distribution in a 
lossless MUT with very high dielectric constant (e, = 20) located on the axis 
of the cylindrical cavity. In Fig. 3.2 the dependencies of the frequency shift Af 
and Q- factor versus e’ and tan 6 are plotted for TMo19 and TE,;, mode in the 
cylindrical resonator. It is obvious that for the TMy;9-mode the frequency shift 
with increasing permittivity shows a nearly linear behavior for 1 < ee, < 50 
without saturation, while for the TE,,,-mode it saturates already at €, > 10. 
Vice versa the quality factor for the TM ;9-mode gets saturated significant 
earlier with increasing tan ô (tan ô > 0.2) as compared to the TE,,,-mode. 
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Figure 3.1: Sideview of cavity E-Field for €i. = 20 and tan 6 = 0.: TE,,, (top left) and 
TMoıo (top right), E-Field stretch in the center of the cavity: TE,,, (bottom left) 
and TMo10 (bottom right) [64] 
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Figure 3.2: 
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Simulated Af(e;.) for tan ô = 0.01 (top) and Q(tan ô) for e, = 40 (bottom) for 
a cylindrical sample with 3 mm in diameter in a cavity with 88.5 mm in diameter 
and a length of 90 mm [64] 


The dual-mode resonator approach solves this problem. The loss factor is 
measured with the TE,,;-mode and the dielectric constant in the TMg10-mode. 
The modes are well separated from each other because of their orthogonality 
(see Fig. 3.1). That allows a separate measurement at both modes without 
interference. Additionally, for samples with a high permittivity the TMo4o- 
mode resonance shows a much larger frequency shift as compared to the TE; 11 
mode resonance. However, because of very low quality factors for tan ô > 0.2 
the TM; 9-mode resonance becomes unmeasurable. This is why the focus is 
on measurement of dielectric losses in the TE,,,-mode resonance. 
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Those effects can be explained with the help of Eq. 2.17 and considering an 
ideal sample without dielectric losses. Fig. 3.1 reveals that for the TE,,,-mode 
the electric field within the sample is decreasing significantly for large €,” 
values. This phenomenon can be explained by the oblique incident electric 
field on the cylindrical MUT surface, which leads to a scattering of the field. 
Therefore, the numerator in the Eq. 2.17 decreases. At the same time, the stored 
energy in the cavity, becomes higher due to the higher effective e; in the cavity. 
Therefore, the denominator becomes higher. This results in a continuous 
decreasing frequency shift and therefore a saturation of the frequency. 


For the TMg10-mode (Fig. 3.1) numerical simulations show that the E-field is 
slightly peaking around the sample position. This can be explained by the 
electric fields oriented along the MUT. Accordingly, the numerator in Eq. 2.17 
is slightly increasing with €, . At the same time, the total stored energy has to 
increase with e,. Then the total stored energy in the cavity outside from the 
sample becomes slightly lower (because of field peaking), that results in 
constant or weak increasing frequency shift function with e;. 


To account the above-mentioned deviations in the dual mode cavity system 
from the classical perturbation theory the calibration of the measurement 
system is realized with the help of interpolated full wave 3D simulation. 
Further information will be given in chapter 3.2.2. These results were 
published in [64]. 


3.1.2 Cavity design 


The designed empty cavity operates at the TE,,,-mode and at the TMy;9-mode 
at around 2.6 GHz [64]. This frequency is chosen in consideration of the 
expected frequency shift (see Fig 3.2) and the operating frequency range of the 
power amplifier. The model of the cavity used for CST simulations is shown 
in Fig. 3.3. The cavity diameter is 88.49 mm and the cavity height is 
88.50 mm. This dimension corresponds to the intersection between the both 
modes, which exist for a diameter equal to the radius (see Fig 2.2). The 
considered dimensions enables the excitation of only the two modes in the 
resonator in the frequency range under consideration. The cavity material is 
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aluminum 6020, which has a specific electrical resistivity of about 3.7 1076 
ohm'em [65]. With the help of Fig. 2.3, the theoretical Q-factor 
(without coupling and under consideration of the cavity material) can be 
determined to be approx. 17500 and 18900 for the TMy,9- mode and the 
TE}11-mode, respectively. The unloaded Q-factor included coupling for the 
TE,,,-mode is calculated to be about 8000 and about 6000 for the 
TMoıo-mode, respectively [64]. 


2a = 88.49 mm 


88.50 mm 


d= 


Figure 3.3: Dual-mode cavity simulated with CST Microwave Studio [64] 


The TE,ı1-mode is coupled through a coupling element at port 1, which is a 
modified N-type coaxial adaptor with a 5.30 mm long copper pin acting as an 
inductive coupling (see Fig. 3.3). The microwave power of the TE,,,-mode is 
decoupled from the cavity at port 2 where a modified SMA coaxial adaptor 
with a 3.05 mm copper pin is used. Through the port 3 the microwave is 
coupled to the TMy;9-mode. The port 3 is a modified N-type coaxial adaptor 
with a copper pin of 11.70 mm in length. The decoupling of the TM); 9-mode 
is realized through port 4, using a modified SMA coaxial adaptor with a 
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3.55 mm long copper pin. The pin length of the coupling ports is selected to 
allow a sufficient strong power coupling for subsequent microwave- assisted 
heating of the MUT, and the decoupling ports are selected to allow a 
measurable transmission (> -60 dB) in the loaded cavity. This coupling design 
enables an isolation of about 60 dB between both modes (see Fig. 3.4). In Fig. 
3.4, the cross-talk between port 1 and port 3 (from Fig. 3.3) is shown, as these 
are the two with the strongest coupling factor. At both resonance positions two 
characteristic ‘fano-like’ resonances are visible [66]. Such resonance results 
from the interference between a broadband, e.g. through the direct transmission 
between the two coupling ports, and a narrow band resonance, e.g. TMg19- 
mode or TE,,,-mode of the resonator. As sample holder a fused quartz-tube, 
with an inner diameter of 2 mm and an outer diameter of 6 mm, is used (see 
Fig. 3.3). The diameter is chosen according to the expected maximum 
frequency shift for a sample permittivity of up to €, = 80 with regards to the 
frequency bandwidth of the amplifier (2.0 — 2.7 GHz). The liquid sample fills 
the fused quartz tube. 


S,, [dB] 


2.2 2.4 2.6 2.8 
Frequency [GHz] 


Figure 3.4: Simulated transmission between port | and port 3 of the dual-mode cavity 
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As already mentioned, here the TMy;9-mode is also used to heat the sample by 
microwaves. This mode is particularly suitable for heating, as it has a 
homogeneous field along the sample. The resulting inhomogeneity in 
temperature along the sample was investigated both in experiment with the 
help of two fiber-optic sensors of the type Optocon TS2 and in simulation as 
well. The result ofthe simulation with COMSOL multiphysics showed that the 
expected temperature inhomogeneity along the sample inside the cavity due to 
convection and heat conduction losses is + 1 °C at around 180 °C (see Fig. 
3.5). This behavior is also verified experimentally. The gradient at both MUT 
ends is due to the transition from the cold MUT part outside the microwave 
field to the heated MUT part inside the microwave field. 


50 100 150 200 
TPC] 


Figure 3.5: Temperature distribution along the MUT by heating with the TMy;) mode and 
taking into account the convection and conduction heat losses calculated with 
COMSOL multiphysics simulations 


3.1.3 Measurement setup 


The block diagram ofthe measurement system is shown in Fig. 3.6 [64]. A key 
element of the measurement system is the frequency-tunable power amplifier. 
The WMP2350/700/50MK-A from G-Way Microwave, which provides up to 
50 dBm (100 W) output power. It has a gain of up to 50 dB. A vector network 
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analyzer (VNA) Agilent N5224A is utilized both for the measurement of the 
transmission coefficient of the TMy10-mode and for providing the input power 
for the amplifier. To reduce the nonlinear effects of elements between VNA 
output and input of the cavity, the so-called “high power” configuration of the 
VNA is used including the external reference signal. A second VNA HP8720d 
is used for measurement of the transmission coefficient of the TE,,,-mode. 
The sample temperature is measured in the center of the sample with the help 
of a fiber-optic temperature sensor of type Optocon TS-NANO. The 
temperature sensor signal is used for the temperature control along the 
predefined temperature scenario with the help of a proportional-integral- 
derivative (PID) control algorithm. The fiber-optic sensor was calibrated 
before each measurement by a one-point calibration at 180 °C with the help of 
a calibration oven from Eurolec of type MK173 and a calibrated thermometer, 
from Mawi-Therm of type P755. 


VNA 
N5224a 


Ref 
Channel 


Cavity 
TE 444 & TMg4o 


VNA 
HP8720d 


Figure 3.6: Block diagram of the dielectric measurements system [64] 
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Fig. 3.7 shows the measurement cell of the developed measurement setup. 
Since the reaction mixture consists of sodium hydroxide solution and PET 
particles, the sample should be oriented horizontally (to avoid sedimentation 
of the solid parts of the sample). The sample holder was fixed and closed at 
both ends with Swagelok connectors, which can withstand a pressure up to 
18 bars. For active control of the cavity temperature two QuickCool Peltier 
elements of type QC-241-1.0-3.9M were installed. This enables the 
temperature stabilization of the cavity to avoid any dilatation, that may 
influence the reference resonance frequency of the unperturbed cavity. 


Figure 3.7: Test setup for the dielectric measurements 


3.1.4 Calibration method 


Following the classical cavity perturbation method, for finding of the dielectric 
constant and loss factor from the experimentally measured frequency shift 
(fo — fi) and (1/Q,; — 1/Qo) the Eq. 2.21 and Eq. 2.22 have to be solved. 
However, for high permittivity values this method is no longer valid while the 
calibration factors A and B are not constant but material-dependent. To resolve 
this problem a novel calibration procedure is used [64]. For that purpose, the 
cavity geometry (see Fig. 3.3) is accurately simulated with the help of 3D full 
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wave simulator CST Microwave Studio (Fig. 3.8). The frequency domain 
solver is used for the range ofthe expected dielectric properties. The frequency 
domain solver is especially suitable for geometries with a very high Q-factor. 
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Figure 3.8: Simulated and measured S,,- parameters for the TE111-mode (top) and for the 
TMo10-mode (bottom) [64] 
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Fig. 3.8 shows the measured and simulated transmission for the two modes in 
the empty cavity. In order to exclude errors due to manufacturing tolerances, 
the exact diameter of the cavity is determined based on the TM,10-mode 
resonance frequency and then the exact length of the manufactured cavity is 
determined based on the TE,,,-mode resonance frequency. 


Then, based on the precisely determined dimensions of the manufactured 
cavity, the perturbed resonator is simulated for different combinations of 
dielectric properties of the MUT. Thus, the Q-factor and f, for the different 
MUT parameters can be determined for both modes. To estimate e’, the 
calculated data of the frequency shift for the TM); 9-mode are interpolated. 
Here a linear interpolation is used, because the frequency shift shows an almost 
linear behavior versus e’. To estimate the loss factor of the TE,,,-mode, the 
change in the Q-factor data is interpolated with a cubic spline. For the 
estimation of tan 6 the a-priori knowledge of €’ is required, since the quality 
factor of the TE,,,-mode depends on e’ strongly. 


3.1.5 Error estimation 


To validate the designed dual-mode experimental setup, two known high lossy 
liquids, alcohol and distilled water, are measured at room temperature. The 
same materials are also measured with the commercial coaxial probe 
HP85070B (Fig. 3.18). The properties measured with both dielectric 
measurements systems and the corresponding literature values are compared 
in Table 3.1 


Table 3.1: Measured dielectric properties with the dual-mode approach, coaxial probe 
HP85070B and the coresponding literature values shown in [43] and [67] 


Material Duel mode cavaly! Coaxial probe Literature [67]; [43] 


Alcohol 10.55 9.75 10.67 9.90 11.23 11.02 
Dist. 
water 
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For both highly lossy liquids, the dual mode cavity approach shows a good 
agreement with the results from the coaxial probe. The literature values are for 
water at 3 GHz [43] and for spirit at 5.65 GHz [67] that may explain the 
deviation from the measured data. 


The most important sources of error are the inaccuracy ofthe sample volume, 
the error in the numerical calibration of the coefficients A and B, the variation 
of the resonator volume and the frequency resolution of the VNA. In the 
following section, the errors are determined individually and the impact on the 
measurement accuracy is estimated. 


Sample volume 


The accuracy of the cavity perturbation method is strongly affected by the 
accuracy of the sample geometry. Since the sample length is equal to the 
resonator length, only the sample diameter can affect the measurement 
accuracy. To minimize this source of error all applied sample holders had been 
cut from the same quartz tube. The inner diameter of the used quartz tube is 
2.05 mm and the maximum measured deviation is + 0.01 mm. This results in 
a measurement error of < 1% for the dielectric constant and < 2 % for the 
loss factor. 


Numerical calibration 


Another uncertainty originates from the numerical calibration of the 
coefficients A and B (Eq. 2.21 and 2.22). The error is caused by the error in 
tolerances of the simulation and the error in the linear interpolation of the 
simulation results. The error due to the interpolation is estimated directly by 
MatLab and the error due to numerical simulations was determined by varying 
the discretization parameters and the pre-set simulation accuracy. The 
maximum error for the coefficients A and B within the calibration ranges 
5 < € < 35 and 0.1 < tan ô < 1.2, amount to 1 % and 2 %, respectively. 


Resonator volume 


Another uncertainty may result from the thermal expansion of the cavity. This 
effect, however, is negligible as the cavity temperature is stabilized to 
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15 +0.3 °C cavity. The cavity was stabilized at 15 °C instead of room 
temperature to ensure a more accurate temperature stabilization, since the 
Peltier elements used can provide an active cooling only. 


VNA resolution 


The limited frequency resolution of the VNA may cause an error. For the 
present measurement system, a relatively low frequency resolution of the VNA 
is used to enable a fast system control. The frequency resolution of the ‘Agilent 
N5224a’, which is used for the measurement of the TMy;9—mode, is set to 601 
frequency points. During measurements, the frequency span is typically 50 
MHz. This results in a maximum error in the measured resonance frequency 
of 50 MHz / 601 ~ 83 kHz. It results in an error of + 0.01 for the dielectric 
constant. In addition, the frequency resolution of the ‘HP8720d’, which is 
utilized for the measurement of the TE,,,; —mode, is set to 601 measurements 
for a frequency bandwidth of 3 MHz. Since the unloaded cavity has a quality 
factor of approx. 8000, it corresponds to a maximum deviation in the Q-factor 
of approx. +60. This corresponds to a maximum error of + 4% for the loss 
factor in the range 0.1 < tan ô < 1.2. 


Estimation of the total error 


In order to define the accuracy of the proposed measurement setup various 
error sources have been investigated independently, as discussed above. 
Assuming that the individual errors are not correlated, the total error can be 
calculated with following equation [68]: 


dy \* 02 dy \* x» 
oy re A 6-1) 
where y is a function of A, ‚..., An. A1» ...,A„ are n uncorrelated variables. 


61; --- O4n are the variances of the corresponding random variables. 


From the Eq. 3.1 and the Eq. 2.21 the error of €; can be estimated to be: 
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The overall accuracy of €; is estimated to be +1.5 %. Similarly, it results from 
equations 2.22 and 3.1: 
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The overall error in the measured tan ô is additionally dependent on the 
inaccuracy in the measured e, , which is prerequired for the estimation of 
losses. From Eq. 3.3 it results in an accuracy of +8 %. 


3.2 Microwave calorimetry 


In the present work, the approach from Nesbitt et al. [37] has been extended 
with advanced numerical simulations to enable a more qualitative thermal 
analysis of the MUT. Through the addition of the thermal analysis, for the first 
time the proposed method enabled the estimation of the specific heat capacity 
of MUT under microwave heating conditions. This system is proposed in [69]. 


The calorimetric analysis bases on the evaluation of the heat flow into or out 
from the MUT during a constant heating rate. In case of microwave calorimetry 
the microwave energy, which remains in the cavity (Paps) is calculated from 
the microwave power balance in the cavity: 


Pabs = Pinc ~ Prefi ~ Prrans (3.4) 


where Pinc and Prerı are the measured incident and reflected power. The 
transmitted power (P;rans) is assumed to be negligible. Furthermore, in order 
to specify the part of Paps, responsible for the measured sample temperature 
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increase, the power measurements is extended by the simulative estimation of 
the wall losses, convection and conduction thermal losses. 


3.2.1 Measurement setup 


Two Anritsu power meters ML2487B and ML2495A are used for the 
measurement of the incident and the reflected power. The isolation between 
the incident and the reflected power is realized by using a bidirectional coupler 
of type Narda-ATM PNR CHP274-30F-30R. For accurate power 
measurements, the VNA and both power meters need a perfect 
synchronization, which is implemented by using a trigger signal generated by 
the VNA. 


VNA 
N5224a 


Trig Out 
WMP2350/700/50MK-A 
Port 2 


Figure 3.9: Block diagram of the power measurements system [69] 


3.2.2 Calibration method 


Calibration of the power measurements 


A power calibration is required in order to obtain accurate values of the 
absorbed power based on Eq. Fehler! Textmarke nicht definiert.3.4. The 
VNA is used to measure the four-port scattering matrix (S-parameters). Those 
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four ports are the ports for the power meters (port 3 & 4), the power amplifier 
output (port 1) and the cavity input (port 2) (Fig. 3.9). The absolute values for 
the incident and the reflected power are calculated using the following 
equations [69]: 


DS, — Bm 


refl 7 
j Sho "Så = S$ Sh 


(3.5) 


.c2 

Ping = a E oz, (3.6) 

31 
P3 and P, are the power levels at port 3 and port 4 respectively measured by 
the power meter (see Fig. 3.9) to monitor the incident and reflected power, 
respectively. The power meter is operated in the pulse-modulated mode with a 
gate time equal to the sweeping time of the VNA. This enables the input and 
reflected power to be determined very accurately by integration over a 


sweeping period of the VNA. 


The measured frequency dependent S-parameters are shown in Fig. 3.10 and 
Fig. 3.11. It should be noted that the isolation between the incident and the 
reflected power is higher than 20 dB. 


2.1 2.2 2.3 2.4 2.5 2.6 2.7 
Frequency [GHz] 


Figure 3.10: S-Parameters with excitation at port 1 of the power measurement system 
from Fig. 3.9 
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Figure 3.11: S-parameters with excitation at port 2 of the power measurement system 
from Fig. 3.9 


In order to specify the part of Paps, responsible for the measured increase in 
sample temperature, full-wave 3D simulations with CST Microwave are done. 
Based on these simulations the cavity wall losses can be estimated accurately. 
Additionally, the exact power distribution in the sample can be approximated. 
For a further improvement of the accuracy of the calorimetric measurements, 
the heat transfer losses from the sample are simulated with COMSOL 
Multiphysics. For the COMSOL modelling, the simulated power distribution 
from the model used in CST MWS is applied as heat source. 


Estimation of the convection and conduction thermal losses 


In this section, an analytical estimation of the convection and conduction of 
the thermal losses of the sample is done for a given sample temperature Ty. 
That estimation is valid for a system which consists of a glass tube with the 
inner radius r; and the outer radius r,. The length of the tube is l (Fig 3.12). It 
is assumed that the tube is completely filled with a liquid and placed in a 
microwave cavity, where a uniform electric field exists. Because of the 
homogeneity of the MUT and of the electric field, the power density is uniform 
inside the tube. 
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resonator 
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Figure 3.12: Resonator model with the correspondings temperature variables 


The heat transfer by heat conduction from the MUT to the outer surface 
of the glass. The total conductive heat flow is given by the following equation 
[70]: 


T-T, 
conductive = aE) (3.7) 
l 


k is the thermal conductivity of the glass tube. T is the temperature of the outer 
surface of the glass tube. 


Heat transfer by convection occurs from the wall of the glass tube to the cold 
air inside the cavity. The heat flux density q” can be calculated if the 
convection coefficient h is known [71]: 


q” = A(T — To) (3.8) 


At this stage, the value of this coefficient is unknown, but from the literature 
[71], it is known that for natural convection cooling it is in the range of 
5 — 20 W/(Km?). 
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The total heat transfer rate is obtained by integrating the heat flux at the surface 
of the glass tube with an outer radius r;: 


2m cl 
q= | | gq’'r.dbdz = 2nr,hl(T — T») (3.9) 
0 0 


The Eq. 3.9 gives the total heat transfer rate in watt transferred from the hot 
tube to the air by natural convection, knowing the temperature at the surface 
of the glass tube T (Eq. 3.10), and the convection transfer coefficient h. 


Assuming that no energy remains stored in the glass tube, the conductive heat 
flow from Eq. 3.7 should be equal to the convective heat flow from Eq. 3.9. 
This allows the temperature of the outer surface of the glass tube T to be 
determined as follows: 


kT, + hr, In (=) To 


T= (3.10) 


hr, In (5) +k 


Tæ corresponds to the temperature of the surrounding air in the cavity and is 
assumed to be equal to the cavity wall temperature. 


The uniform absorbed power Pyuyr by the MUT is given as follow: 
Puur = mc,AT/At + 2nr,hl(T — To) (3.11) 


m corresponds to the mass of the MUT. c, is the specific heat capacity of the 
MUT. The first term from Eq. 3.11 represents the absorbed power needed to 
elevate the material temperature by AT in time At, the last term is the power 
transmitted to the surrounding air inside of the cavity by convection. From Eq. 
3.11, the temperature rise AT can be calculated: 


mc 


AT (3.12) 


p 
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COMSOL model 


Using a 2D axisymmetric model, a coupled fluid and thermal model is solved 
to obtain the convection flux and the transfer coefficient h. The model solves 
the laminar Navier-Stokes equation [71] together with the heat equation. The 
convection takes place due to the differences in the air density p, caused by 
the temperature gradient. The volume force f(r,z,t) in COMSOL is 
implemented as follows: 


fC, zt) = g` (Pa(7,Z,t) — Pref) (3.13) 


r and z represent the cylindrical coordinates and t is the time. 
Pref = 1.19 kg/m’, is the density of air at an ambient temperature of 293 K 
and at a pressure of one bar. The volume force f at the location (r, z) and at the 
time t is proportional to the current local density. The earth gravity is given as 
g. The air is modelled as an ideal gas, and the density is expressed as a function 
of the temperature and pressure via the ideal gas law [72]: 


Ma R 
PaVa = NaRT = mM eT = Ma T =MaRal . (3.14) 
a a 
Ma 
Pa = y Rat = PaRaT (3.15) 
a 
Pa 
Spa RT (3.16) 


R = 8.3 J/(mol-K) is the universal gas constant and M, = 29 g/mol is the 
molar mass of dry air. The specific gas constant for air is assumed to be 
Ra” 286 J/(kg K). 

Geometry 


The inner and outer radius of the glass tube are r; ~ 1.5 mm and r, = 3 mm 
respectively. The diameter and height of the cavity are 88 mm each. 
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Material properties 


For the glass tube is assumed a density ofp = 2210 kg/m? [73]. The specific 
heat capacity c, of the quartz tube is assumed to be 770 J/(kg K) and the 
thermal conductivity is assumed to be 1.4 W/(m K) [73]. The thermal 
properties of air are a function of temperature and follows the ideal gas law. 


Power Calibration 


The power calibration is performed in the following manner: firstly, the 
dielectric properties of the MUT and the incident and the reflected power are 
measured. Secondly, based on the measured dielectric properties, the electric 
field density and the wall losses in the cavity with CST microwave studio are 
estimated. Thirdly, the simulated data are used as input for the multi-physic 
modelling using COMSOL (see Fig. 3.13). Based on the material and geometry 
parameters mentioned above, the simulative estimated transfer coefficient h is 
13 W/(m? K). 


Measurement 
(er, tan(ô), Pinc» Pref) 


| Comsol multiphysics: heat losses 
through convection & conduction 
CST MWS: electric field u 
a 
Results 
Pabs,ert(t T) 


Figure 3.13: Power calibration methodology 
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3.2.3 Validation of the system performance 


To validate the measurement system three known liquids have been measured. 
A sample with a radius of 1.5 mm is used to keep the power requirement as 
high as possible. In this way the absolute error on the heat capacity 
determination can be kept low since the energy required to heat the sample is 
much greater than the thermal or wall losses. All measurements of the specific 
heat capacities are made at a constant power level of 25 Watt for the incident 
power and a fixed frequency span of 0.1 MHz around the resonance frequency. 
In this way the influence of temperature control is minimized while 
maximizing the absorbed power, thus minimizing the relative error of the 
power measurements. Only in this way a determination of the heat capacity 
was possible. 


Using water as a first sample, this results in an average absorbed power of 
5.71 Watt. According to the simulations the wall losses amount to 
approximately 31 % of the absorbed power. From Eq. 3.9 the average 
thermal losses are calculated to be approx. 0.5 Watt. This results in an absorbed 
power of about 3.44 Watt, which is responsible for the heating of the water 
sample. From Eq. 3.11 the measured specific heat capacity is approx. 5.5 J/(g 
K). 


The measured properties from the three measured liquids with the developed 
calorimeter and the literature values [74] are presented in Table 3.2. 


Table 3.2: Comparison between measured and literature values of the specific heat capacity 


Material Cp [74] Cp (measured) Cp (corrected) 
[J/(g:K)] [J/(g:K)] [J/(gK)] 
Ethanol 2.4 3.2+0.5 2.44 
Ethylene glycol 2.2 2.8 + 0.4 2.14 
Water 4.2 5.5+0.7 4.2 


For all three liquids, the proposed microwave calorimetry shows the right order 
of the specific heat capacity. The most important sources of error are the 
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inaccuracy of the measured dielectric properties, the inaccuracy of the sample 
and cavity temperature measurements, the measured incident and reflected 
power and finally the error in the numerical calibration. 


Table 3.2 shows a systematic overestimation of the measured heat capacities. 
A further improvement of the heat capacity determination can be achieved by 
an extended calibration of the measurements system by a known sample 
(e.g. water). This correction can be made by multiplying the measured specific 
heat capacity by the ratio between the actual and the measured value of the 
well-known material (reference material). In case of Table 3.2 and assuming 
that the heat capacity of water is well known and equal to 4.2 J/(g:K), resulting 
in a correction factor of ~0.76, the measured heat capacity of ethanol and 
ethylene glycol can be corrected to 2.44 J/(g-K) and 2.14 J/(g-K), respectively. 


An important goal of calorimetric measurements is to find out the exothermic 
or endothermic behavior of the depolymerization reaction. Hence, it has to be 
investigated for known reactions, e.g. for the curing of epoxy resins. For this, 
a quarz tube with an inner diameter of 3 mm is completly filled with the epoxy 
resin ‘Pattex Stabilit Express’ [75]. On the basis of a material density of 
1.172 g/cm? this amounts to a material quantity of 0.7357 g assuming a 
volume of 0.627 cm?. Fig. 3.14 shows the curing process of the epoxy resin 
with two different constant heating rates. More exactly, it shows the absorbed 
power measured during the curing process of an epoxy resin. An approximate 
linear increase of the absorbed power demand with increasing temperature can 
be recognized, which is due to an increased heat loss from the sample surface 
to the cold environment by conduction and convection. As soon as the 
exothermic polymerization reaction starts at about 80 °C, the power 
requirement to maintain a constant heating rate decreases until the 
polymerization is completed. Furthermore, as the heating rate increases, the 
observed decrease in power demand during polymerization increases. In ideal 
case, the Fig. 3.14 provides a proper estimation of the energy released during 
the reaction. That energy can be calculated by integrating the surfaces ‘area 1’ 
and ‘area 2’. The calculated energy is 416 J and 392 J from the ‘area 1’ and 
‘area 2’, correspondingly. This corresponds to a released energy during the 
exothermic curing reaction between 532 J/g and 565 J/g, which is comparable 
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to literature values of about 400 J/g for similar epoxy resins [76]. Taking water 
as reference material, the measured released energy during the reaction can be 
corrected to 406 J/g and 432 J/g. 
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Figure 3.14:  Calorimetric monitoring of curing of epoxid resins 


The power measurements should fulfill two main goals. First, it should make 
an exothermic or endothermic behavior visible and second, the heat capacity 
of the investigated material should be determined. The inaccuracy of the power 
measurements should be investigated for systematic and random errors. 
Hereby, the systematic errors only affect the determination of the heat 
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capacity, since relative measurements are sufficient for the determination of 
the behavior of the reaction. The main error sources in the estimation of the 
incident and the reflected power can be calculated with the help of the 
equations 3.5 and 3.6. It follows that the errors are caused by the inaccuracy of 
the measured power signals P} and P, with the power meter and the inaccuracy 
ofthe measured S- parameters. According to [77] the systematic errors of the 
power sensors are about +4 % and the random errors are +0.11 %. The error 
of the measured S-parameters is about +0.01 dB [78]. It can be considered as 
a purely systematic error. These errors can be calibrated using a well-known 
sample (e.g. water). From Fig. 3.10 and 3.11 it can be seen that the transmission 
coefficients S32 and Sy, are significantly smaller than the transmission 
coefficients S3, and S42. This simplifies the equations 3.5 and 3.6 as follows: 


Pe (3.17) 
ref Ss 
S. 2 
Pine © (=) Pz (3.18) 
31 
The equations 3.1, 3.17 and 3.18 result in 
op, ôS,» 
Fa 1 
Prefl P? +4 5, G 9) 
65, ö2 62 
a +42 +42 (3.20) 
RE fe Sh S31 


For the calculated Pinc and Prerı this results in a systematic error of about 
+4.1 % and +4 %, respectively. Since the power levels of Pinc and Prerı are 
about 10 Watt, an error for Paps of up to 0.8 Watt is feasible and accordingly 
the determination of the heat capacity is very imprecise. The random error in 
the calculation of Pinc and P,erı is about +0.11 % and a random error for Paps 
of up to 0.022 Watt is feasible. It follows that the power measurements can be 
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considered as meaningful in the determination of the endothermic/ 
exothermic behavior of the reaction. 


3.3 Combined dielectric and calorimetric 
measurements 


3.3.1 Measurements setup 


Finally, both subsystems for the dielectric and calorimetric measurements were 
merged. This enables the simultaneous temperature dependent dielectric and 
calorimetric characterization of materials under microwave heating conditions. 
The block diagram of the measurement system is shown in Fig. 3.15. 


VNA 
N5224a 


Trig Out 


Meter 


ML2495A | 


VNA 
HP8720d 


Figure 3.15: Block diagram of the measurements system 


3.3.2 Automatic calibration and control system 


Operation control and data acquisition are performed on a personal computer 
(PC) with the help of a MATLAB code. In the code, the standard commands 
for programmable instruments (SCPI) are used to control the I/O interfaces. 
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The network analyzer HP9720d is controlled via a remote general-purpose 
interface bus (remote GPIB). The second network analyzer N5224a is 
controlled via a local area network (LAN), as well as the two power meters. 
The power supply HCS-3202 control and the acquisition of the thermocouple 
data (NI9211) are controlled via USB ports. The Optocon fiber optic 
temperature sensor is controlled via RS-232 port. The software allows a loop- 
controlled operation including heating and data acquisition (see Fig. 3.16). 


Initialization 
Cavity at 15 °C 


Read Pinc and Prefi 
Center peak S,} at 
HP8720d 
Read f Qat 
HP8720d 


Calculate €,,tanð 
and save data 


Read MUT 
temperature 
PID controller 

Center peak S31 set 
span at N5224a 


Read f, Q at N5224a 


Cavity at 15 °C 


Meas. 
completed? 


Meas. 
completed? 


Figure 3.16: Scheme of the control of the dielectric and calorimetric measurements set up 
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The achievable control time is 0.5 s. In this setup the main restricting factor for 
the control time is the communication with the VNA, which is about 0.45 s. 
For this reason, the main code is parallelized. It consists of two well 
synchronized loops. The temperature stabilization of the resonator works 
independently during the measurements process. The structure of the 
developed control algorithm is shown in Fig. 3.16. Furthermore, a graphical 
user interface (GUI) for the implementation of fully-automated measurements 
has been developed. 


The temperature control is carried out with the help of a proportional-integral- 
derivative (PID) controller. The PID controller uses as input the measured and 
the set temperature values and provides a frequency step as the output, which 
is added to the actual frequency. This span value is given as input to the VNA 
N5224a. The smaller the frequency step the larger the heating rate of the 
sample and vice versa. 


The equation for the PID controller can be formulated as follows [79]: 
Aspan ‘= Kp * ((e —@,)+K,:e+Kp:(e-2e, +e, )) (3.21) 


where K,,K;,Kp are the proportional, integral and derivative gains, 
respectively. e, e,,e, are the actual, previous and two times before error value, 
respectively. The error value is defined as the difference between a nominal 
temperature Tyominaı and an actual temperature Tactuaı- For the developed 
measurement setup the following control constants were determined 
experimentally with the help of a multi-step temperature scenario (see Fig. 
3.17): Kp = 0.65, K; = 0.28, Kp = 0. 
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Figure 3.17: Example of temperature control with optimized PID controller parameters 


The starting parameters of the PID controller are selected in a way that an 
overshooting at the beginning of the process (close to room temperature) is 
avoided. Otherwise, the settling time of the PID controller would be very long 
due to the lack of an active cooling. Fig. 3.17 shows that the higher the actual 
temperature is the shorter the settling time and the lower the overshoot 
becomes. This can be explained by the increasing conduction and convection 
heat losses due to the cold environment. That enables a faster and more 
accurate reaction to the temperature by variation of the microwave power. 


In ideal case, the short control time in combination with the PID controller 
allows accurate temperature control with a maximum deviation from the preset 
temperature profile of about +0.3 °C. But as the dielectric loss tangent 
becomes larger (tan ô > 0.5), the temperature control becomes more and more 
complex. Because of the very high dielectric losses the microwave cavity 
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develops a very low-quality factor and therefore becomes very broadband. 
Accordingly, the variation of the microwave power and the corresponding 
temperature control can be realized in a very restricted way, since all the 
temperature control is done by the frequency span variation. In that case, 
greater temperature fluctuation of up to +3 °C are common. 


3.4 Measurements using the coaxial probe 


In the present work, the various reactants are partially measured depending on 
temperature and frequency by using a commercial coaxial probe of type 
HP85070B [80]. The coaxial probe enables the dielectric characterization of 
materials from DC up to 20 GHz in the temperature range from -40 °C to 
+200 °C. The typical accuracy is +5 % for the dielectric constant and +0.05 
for the dielectric loss tangent [81]. A magnetic hotplate stirrer VWR VMS-C7 
is used for the controlled heating of the liquid probe up to 100 °C. A picture of 
the measurement setup is shown in the Fig. 3.18. 


Thermocouple HP85070B 


Figure 3.18: Measurement setup based on commercial coaxial probe HP85070B 
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4.1 Sample preparation 


For the following measurements, MEG and PTA with a purity grade of 
> 99 % and = 98 %, are obtained from company Alfa Aesar. NaOH with a 
purity of > 99 %. is obtained from company Merck Millipore. The PET 
particle size is about 500 um (see Fig. 5.13). 


Table 4.1: Mixture composition at the different reaction steps 


Reaction 

es 

| enfosi- foantses [1830 140 17.78 | 82.4 | [5241206] 2.06 | 
ii 


30 pelose a67 1038| 97> [1507] 
40 
50__| 90 | 0.47 | 77.8 [0.47| 984 | 15.85 | 
60 
70 __| 54 | 0.28 | 108.9] 0.66| 995 | 16.04 | 
80 
90 __| 18 | 0.09 | 140.0] 0.84 | 1007 | 16.23 | 
100 | 0 | 0.00 | 155.6] 0.94] 1013 | 16.32 | 


In Table 4.1 different chemical mixtures are given. Each mixture represents a 
different level of the reaction progress, e.g. 0 % is equal to a status of no 
depolymerization, whereas 100 % is equal to full depolymerization. As 
expected by eq. 2.23, one mol PET reacts with two mol NaOH. It yields one 
mol PTA and one mol Na,PTA. Since Na,PTA is not directly available as a 
substance, here NaOH and PTA are mixed. It yields to Na,PTA. 
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4.2 Dielectric measurement results 


4.2.1 Dielectric properties of the reaction mixture 


The samples have been prepared just prior to the individual measurements, to 
avoid any reaction before the actual measurement. The measurement 
procedure including constant heating up to 180 °C and subsequent controlled 
cooling did last approximately one hour. Eleven samples are tested with a 
depolymerization level of 0 % up to 100 % as shown in Table 4.1. 


In Fig. 4.1, the time- and temperature-dependent dielectric properties of the 
reaction mixture with PET/PTA are presented. During the experiment, the 
heating rate was chosen to be 3 °C/min and the cooling rate -5 °C/min. Larger 
heating rates were not feasible because of the long reaction time of the 
thermo-sensor, the risk for boiling within the sample volume and leaking 
problems of the sample from the sample holder. During the initial experiments 
such problems were noticed as well as the local heating of the inner glass tube 
surface which was affecting the glass itself in a way that it could not be cleaned 
any more. 


Fig. 4.1 shows the temperature-dependent properties during heating and 
cooling of the samples. The direct comparison of the heating and cooling 
scenarios for 0 % depolymerization (red curve) reveals information about the 
influence of the chemical process onto the dielectric constant and dielectric 
loss. Both show higher values during heating as compared to the cooling 
scheme. The sample ‘100 % depolymerized’ (green line), which represents the 
dielectric properties of finished chemical depolymerization, follows quite well 
the cooling curve of sample ‘0 % depolymerized’, indicating that the chemical 
process is finished as soon as the target temperature of 180 °C is reached. This 
is probably due to the slow heating rate of 3 °C/min as well as by the fact that 
the PET used was a fine powder. In case of faster heating or larger PET 
particles still some isothermal dwell at the targeted reaction temperature might 
be needed. Due to the low heating rate, the direct characterization of the 
mixtures at a process temperature of 180 °C was not possible since the 
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chemical reaction already started during the heating and obviously finished 
before the target temperature was reached. 
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Figure 4.1: Temperature depending permittivity for all samples according to Table 4.1 


Three different effects that take place during heating can explain the main 
course of the dielectric parameters from Fig 4.1. First, the dielectric properties 
of MEG (see Fig 8.1) dominates, which also forms a large volume fraction 
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(~74 %) of the total mixture. This factor is mainly responsible for the strong 
reduction of tan ö at the beginning of the heating phase and the simultaneous 
increase of €,. At the same time, the conductivity losses caused by NaOH 
increases with the temperature (see Fig 8.2). Due to these losses, the slide 
increase of tan 6 can be explained by temperatures T > 80 °C). And finally the 
depolymerization reaction, which leads to an effective reduction of the NaOH 
concentration. The effective reduction of the concentration of unbound NaOH 
may be responsible for the reduction of tan 6. 


For the latter system design the permittivity measured for the sample ‘0 % 
depolymerized’ is most relevant since it represents the dielectric properties 
along the full process. The corresponding properties for all samples according 
to Table 4.1 at room temperature and at 180 °C after heating with a heating 
rate of 3 °C/min are listed in Table 4.2. 


Table 4.2: Measured dielectric properties for the mixtures listed in Table 4.1 


Reaction 

progress €-@ 25°C | &@ 180°C | tan @25°C| tan @ 180°C 
[%] 
0 14.9 21.2 0.95 0.47 
10 14.7 20.8 0.91 0.48 
20 14.4 20.9 0.85 0.48 
30 14.2 20.9 0.82 0.48 
40 13.8 21.0 0.76 0.48 
50 13.5 21.0 0.71 0.48 
60 13.1 21.1 0.67 0.48 
70 12.6 21.1 0.61 0.48 
80 12.0 21.1 0.55 0.48 
90 11.7 21.2 0.51 0.48 
100 11.6 21.1 0.49 0.48 


It must be mentioned, that all samples are 100 % depolymerized as soon as the 
target temperature of 180 °C is reached. This explains why all samples show 
the same permittivity at 180 °C. A future applicator must be able to work 
efficiently over the full range of permittivities, which means for a 
permittivity 15 < €, < 25. The difference between the heating and cooling 
scenario for 100 % depolymerization (green curve) may be explained at least 
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partially by a temperature demand to enable the binding of NaOH and PTA 
to Na$PTA. This can be investigated by the permittivity curve of the re-heated 
sample. A re-heating of the sample is only possible up to approx. 100 °C, 
because the chemicals damage the glass fiber sensor after a certain time and a 
repair/replacement becomes necessary. From Fig. 4.2 is evident that the 
direction of the permittivity during the cooling down process is very similar to 
that of the re-heated reaction mixture. 
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Figure 4.2: Temperature dependence of €;. for sample ‘Reaction progress 100 %’ 
according to Table 4.1 incl. re-heating 


An important parameter that can act as an indicator for a maximum material 
thickness is the penetration depth D,. [82]. Based on the measured dielectric 
properties (Fig .4.1) and the Eq. 2.6 the penetration depth during the reaction 
can be calculated. Fig. 4.3 shows the behavior of the penetration depth with 
temperature and reaction progress. As can be seen in the graph the lowest 
penetration depth is estimated to be about 5.5 mm at 25 °C. 
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Figure 4.3: Dp of the sample “Reaction progress 0 %’ 


4.2.2 Dielectric properties of the initial mixture 


For the modelling of the reaction, the temperature-dependent loss factor of the 
reactant (initial mixture), the reaction mixture (initial mixture during the 
reaction), and the reaction product (final mixture) as well as the time- 
dependent dielectric losses of the reaction mixture have be to be known. 


As mentioned in section 4.2.1, it is not possible to characterize directly the 
temperature-dependent properties of the reactants, because the reactants start 
reacting already at low temperatures. At the same time a very fast heating of 
the mixture is not possible due to the slow temperature sensor. Therefore, the 
temperature-dependent ‘initial mixture’ is measured in two steps, separating 
the reactant from the PET. In this way, a reaction during the heating process 
can be avoided. Firstly, the mixture consisting of mono ethylene glycol, water 
and NaOH is measured up to 180 °C. Secondly, the loss factor from virgin PET 
is measured separately up to 180 °C. 


For the dielectric characterization of PET a measuring system based on the 
TE,,;,—Mode was used. The detailed description of this measurement setup can 
be found in [83]. A cylindrical sample with a diameter of 7.8 mm and a length 
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4.2 Dielectric measurement results 


of 10 mm was used for this measurement. The measured dielectric properties 
are shown in figures 4.4. 
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Figure 4.4: Dielectric properties of PET 


For the characterization of the reaction mixture without PET the dual mode 
system is used. The measured dielectric properties are shown in figures 4.5. 
The highly noisy dielectric parameter of MEG, H>O and NaOH without PET 
in Fig. 4.5 is explained due to the very high dielectric losses (tan 5 ~1). In this 
range, both the temperature control and the measurement accuracy of the 
system reaches its limits. In any case, these results are sufficient for further 
consideration. 
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4 Dielectric and calorimetric measurement results 
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Figure 4.5: Dielectric properties of MEG, H,O and NaOH 


From the figures 4.4 and 4.5 is evident that the permittivity of PET is much 
smaller than that of the rest of the reaction mixture. Accordingly, PET does not 
significantly influence the course of the permittivity of the unreacted mixture. 


The temperature dependent dielectric properties of the mixture are calculated 
using Eq. 2.8 with the LR exponential factor ¢ = 1. The mixing formulas for 
the dielectric properties of the full reaction mixture are as follows: 


€, = Vi-per€ri-per + (1 — Vı-per)Erper (4.1) 
€, = Vi-per€ri-per + (1 — Vı-per)Erper (4.2) 
tan 6, = eyı /erı (4.3) 


V,-pgr is the volume percentage of H20, NaOH and MEG in the mixture at the 
beginning of the reaction. The dielectric constant of the full-unreacted reaction 
mixture, reaction mixture except of PET (consisting of H,O, NaOH and MEG) 
and PET are e/,,E;,_prr and €; pr, respectively, and €;),,€;) per and ErpET 
are the corresponding loss factors. The volume percentage of PET in 
the reaction mixture is 11 %, from a PET density of 1.38 g/cm? [84], MEG 
density of 1.11 g/cm? [85] and NaOH density of 2.13 g/cm? [86]. The 
calculated tan 6 for the full reaction mixture, based on measurement results 


64 


4.2 Dielectric measurement results 


given in Fig 4.4 and 4.5 following eqs. (4.1) to (4.3) is shown in Fig. 4.6. This 
process can be described similar to the process shown in Fig. 4.1, based on the 
dielectric losses of MEG (Fig 8.1) and NaOH (Fig 8.2). 
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Figure 4.6: tan ô of the unreacted mixture consisting of MEG, H20, NaOH and PET 


4.2.3 Dielectric properties of reactants measured with 
coaxial probe 
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Figure 4.7: Measured dielectric loss tangent of MEG for six different temperatures. 
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4 Dielectric and calorimetric measurement results 


Due to the very high dielectric losses of NaOH and MEG, the dielectric 
properties of the liquid components of the reaction product have been 
measured at different temperatures with the coaxial probe HP 85070B (see 
section 3.4). 


Since the reaction mixture consists mainly of MEG, the loss tangent of MEG 
is firstly considered separately for six different temperatures (30 °C, 40 °C, 
50 °C, 60 °C, 70 °C and 80 °C). The measured data is shown in Fig. 4.7. 


Furthermore, due to the expected high conductivity of the NaOH solution [87], 
the NaOH concentration influences the total permittivity strongly. For a more 
detailed analysis, a separate consideration of the NaOH solution has been 
realized with the coaxial probe method. In each measurement series a different 
concentration of NaOH (5 wt.%, 7.5 wt.%, 10 wt.% and 15 wt.%) is used and 
the complex permittivity is measured for the six different temperatures (30 °C, 
40 °C, 50 °C, 60 °C, 70 °C and 80 °C). 


The weight fraction of NaOH in water is defined as the NaOH mass divided 
by the total mass. 


In Fig. 4.8 the dielectric losses versus temperature for the different NaOH 
concentrations at 2.45 GHz are shown. According to Fig. 4.8, the losses are 
related to the temperature. An increase in the concentration of NaOH, leads to 
a significantly greater temperature dependence of the losses and higher losses. 
This phenomenon has already been described in [88]. 
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4.2 Dielectric measurement results 
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Figure 4.8: NaOH solution: loss tangent versus temperature for four NaOH concentrations 


In addition, the dielectric loss tangents of the entire mixture consisting of 
NaOH, H,O and MEG is characterized. Thereby a fixed weight ratio between 
NaOH and water (37 wt.% according to Table 4.1) was used. The 
concentrations mentioned below are to be understood as weight of NaOH in 
relation to the weight of NaOH, H,O and MEG. In Fig. 4.9 the dielectric loss 
tangents versus temperature for the different NaOH concentrations at 2.45 GHz 
are shown. 
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Figure 4.9: NaOH, H,O and MEG solution: loss tangent versus temperature for four 
NaOH concentrations 
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4 Dielectric and calorimetric measurement results 


In order to derive a mathematical expression for describing the concentration 
dependency, a fitting routine is used. This provides an approximation for the 
losses over the NaOH concentration for a selected frequency at room 
temperature (see Eq. 4.4). Here too, the weight proportion between NaOH and 
water was fixed at 37 wt.%. The resulting progression at room temperature for 
different concentrations is shown in Fig. 4.10. From Fig.4.10 it is visible that 
the dielectric losses follow a path that reaches its maximum for a concentration 
of approximately 15 wt.%. Similar behavior of the NaOH conductivity is 
reported in [89]. The dielectric losses across the NaOH concentration can be 
fitted with a polynomial of third order (see [89]). A comparison of the 
approximated and measured parameters in Fig. 5.11 shows a very good match. 
The expected error through these fitting is less than 2 %. 


= 3 2 
tan Ö,_per,rr = 43 CNaon + A2Cnaon + AıCnaou + Qo (4.4) 


where the coefficients a, = 5.425 -10”°,a, = —4.392 : 1073, a} = 0.1012 
and a, = 0.5912 
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Figure 4.10: Measured and fitted tan ô from NaOH with water and MEG at room temperature 


From Eq. 2.23 it becomes evident that during the alkaline hydrolysis of PET 
one mole of PET reacts with two mole of NaOH. During the reaction NaOH 
binds to Naz[PTA]. This effect leads to an effective reduction of the NaOH 
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4.3 Calorimetric measurement results 


concentration. The dielectric losses of NaOH are disproportionally higher than 
those from the total mixture. This is why a prevalence of the dielectric losses 
throughout the reaction from NaOH is expected. 


In the literature, the molar weight of PET and NaOH are given as 192.2 g/mol 
and 39.997 g/mol, respectively. Under the assumption that the effective 
dielectric loss during the reaction changes due to the precipitation of disodium 
terephthalate and resulting reduction of ionic conductivity, the loss tangent 
development can be modelled based on the data in section 4.2. From Eq. 2.23 
it is obvious that the degree of conversion is linearly proportional to the change 
of NaOH concentration. Depending on the reaction stage the NaOH 
concentration Cyaoy changes between 6.07 and 0.55 wt.% of the total reaction 
mixture. Therefore, the relationship between the degree of reaction and the 
NaOH concentration can be described as follows 


Cnaon = 5.52 (1 — æ) + 0.55 (4.5) 


4.3 Calorimetric measurement results 


In Fig. 4.11, the time-dependent absorbed power of the reaction mixture with 
PET/PTA is plotted. The heating rate is 3 °C/min. Those measurements are 
realized in parallel to the measurements from section 4.2. Due to the slow 
heating rates and the permittivities measured in section 4.2, it can be assumed 
that the reaction mixture is completely reacted during the heating phase. 
Hence, only the heating phase is considered for the further analysis. 


The absorbed power shows an irregular progression for temperatures 
< 30 °C. This can be explained by the time the controller needs to adjust 
the heating parameters. The direct comparison of the heating scenarios for 
‘0 % depolymerization’ (blue curve) and ‘100 % depolymerization’ (orange 
line) is expected to reveal information about the exothermic and endothermic 
behavior of the chemical processes, respectively. The MUT weighted 
approximately 0.335 g. The integration of the gray area ‘Area 1’ in Fig. 4.11 
results in released energy from ~526 J, which under correction with the 
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4 Dielectric and calorimetric measurement results 


reference material water to a released energy of about -368 J follows, that 
corresponds to released energy of ~1591 kJ/mol. This value is comparable 
with the values referred in [58]. Due to the uncertainty of the power 
measurements, further modeling of the reaction kinetics is realized based on 
the permittivity measurements. 
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Figure 4.11: Time dependence of absorbed power during the depolymerization of PET and of 
the reaction product according to Table 4.1 
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5 Modelling of the reaction 
kinetics based on tan 6 


This chapter describes a novel modelling approach to describe the loss factor 
of the reaction mixture during the alkaline hydrolysis of PET, which takes into 
account a progressive decrease in the surface area of the PET particles, in 
contrast to comparable studies [90]. 


5.1 Data preparation for modelling 


For the modelling of the reaction kinetics, the data of the reaction product and 
of the reaction mixture during the reaction are needed. The measured 
tan 6 data of the reaction product (Sample ‘Reaction progress 100 %’) versus 
temperature can be seen in Fig. 5.1. The corresponding fitted data used for 
the further modelling can be found in the Appendix. 
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Figure 5.1: Loss tangent of the reaction product 
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5 Modelling of the reaction kinetics based on tan 6 


Fig. 5.2 and Fig. 5.3 show measured data in dependence ofthe temperature and 
time during the reaction under a heating rate of 3 °C/min, respectively. These 
data correspond to the sample ‘Reaction progress 0 %’ from Table 4.1. Fig. 5.3 
shows a plateau in the time range < 400 s (red circle). This is caused by the 
settling time ofthe PID controller and has nothing to do with the reaction itself. 


1 T 
— Measured 


“20 40 60 80 100 120 140 160 180 
TPC] 


Figure 5.2: Temperature dependence of the loss factor during the reaction (dT/dt = 3 °C/min) 
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Figure 5.3: Time dependence of the loss factor during the reaction (dT/dt = 3 °C/min) 
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5.2 Reaction model 


In order to enable a calculation with Matlab, the measured data were described 
for later modelling by analytical functions. Fig. 5.4 shows the fitted curves of 
the mixtures at their different stages: unreacted (solid black line), reacted 
(dashed black line) and during the reaction (dotted black line) (see 
Appendix A). The red area corresponds to the portion of the mixture that has 
reacted and the green area to the portion ofthe unreacted mixture. 


— -reaction product 
— unreacted mixture 
| a reaction 
Wi rcacted 
0.8 - El unreacted 
0.7 
wo 
c 
S 0.6 
0.57 
0.47 
0.3 
0.2 
20 40 60 80 100 120 140 160 180 
TC] 


Figure 5.4: Temperature dependence of the tan ô of the unreacted mixture, reaction product 
and during the reaction 


5.2 Reaction model 


By substituting Eq. 4.6 in Eq. 4.5, the tan ö,_pgr,rr at room temperature can 
be calculated as follows 


tan ÔL-PET,RT 


= a; (5.52(1 — æ) + 0.55)? + a,(5.52(1 — a) + 0.55)? (5.1) 
+ a,(5.52(1 — a) + 0.55) + ay 
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5 Modelling of the reaction kinetics based on tan 6 


For the modelling, Eq. 2.9 is used to mix tand,_pgr,rr calculated using 
Eq. 5.1 with tan öpgr (see Fig. 4.4). The fractional volume of PET Vpgr is 
assumed to be ~11 % (see Table 4.1). It is also assumed that the density and 
dielectric losses of PET and PTA are the same. From Eq. 2.9 the total loss 
tangent at room temperature can be calculated as follows 


tan ôL rr = Vi-per * tan ö,_per,rr + (1 — VL-per) * tan Oper rr (5.2) 


In Fig. 5.5, the modelled development of tan ô versus the degree of the reaction 
(with the help of Eq. 5.3) is presented. A good agreement between the modelled 
and measured parameter (see Table 4.2) is evident. The maximum deviation 
between the measured and modelled tan ö is < 3 %, which leads according to 
Table 4.2 to an incorrect estimation of the actual reaction stage of up to 6 %. 


1 T T 
* e measured 
— modelled 


10 20 30 40 50 60 70 80 90 100 
a [%] 


Figure 5.5: Measured and simulated tan ô; gr during the reaction at room temperature 


Under the assumption that tan ô behaves similar at room temperature (see 
Fig. 4.8), tan ô of the reaction mixture can be easily estimated over T with the 
help of Eqs. 5.2, 8.3 and 8.4. 
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5.2 Reaction model 


tan 6, (5.3) 
= tan Ornd + (tan ÔL RT — tan Ornarr) 


tan Östart — tan Sena 


tan Östart,rr — tan Ögna,RT 


tan Östare,rr and tan Özna,rr are the dielectric losses of the unreacted reaction 
mixture and the reaction product, respectively. These can be easily calculated 
from Eqs. 8.3 and 8.4 for T = 20 °C. In this way, the measured tan 6 can be 
converted to the corresponding reaction degree. 


The corresponding temperature dependent tan 6, of the reaction mixture over 
the reaction progress is shown in Fig. 5.6. 
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Figure 5.6: Temperature dependent tan 6, of the reaction mixture over the reaction 


The reaction degree a results from a continuous integration over time. In 
general, the reaction degree can be derived from the following differential 
equation [91]: 
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5 Modelling of the reaction kinetics based on tan 6 


Aa 
er 5.4 
me er 
f as function of a comprises the reaction model and x is the temperature 
dependent reaction rate that is given by the Arrhenius formula [92]. 


Kadem) (5.5) 


R = 8.3144598 J/(mol K) is the universal gas constant. A is the Arrhenius 
constant and E, is the activation energy. The latter is the energy required to 
start a chemical reaction. 


Depending on the type of reaction or the function of the reaction mechanism, 
f(«) can be expressed by different formulas [91]. In case of the alkaline 
hydrolysis of PET, f (æ) has been reported as follows [90]: 


f = Ns V Pr Nor" (M — 2a)Pr (5.6) 
where 


e q represents the degree of the PET conversion, 
e Nasis the surface area of the PET flakes, 

e V is the total reaction volume, 

e Nj; is the moles of PET repeating units, 

e M is the molar ratio of NaOH and PET, 

e and, finally b, is the reaction rate order. 


The order of the reaction rate refers to the power dependence of the rate on the 
concentration of the individual reactants [93]. In a first-order reaction, the rate 
is therefore dependent on the concentration of the single species. A second- 
order reaction refers to a reaction whose rate depends on the square of the 
concentration of a single reactant or the combined first-order dependence of 
the rate on the concentrations of two different reactants [93]. According to [90], 
in case of alkaline hydrolysis of PET the reaction rate b, is a first order 
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5.2 Reaction model 


reaction. Here, for the further calculations b, is set equal one. Accordingly, 
from equations 5.4, 5.5 and 5.6 follows. 


= - K+ Nas * (M — 2a)/V (5.7) 


In Eq. 5.7 all parameters are known or measured except K. 


The PET mass during the reaction changes proportionally to 1—a@. The 
volume of PET Vpgr is calculated from the mass divided by the density of 
1.38 g/cm? [84].For further calculation, it is assumed that the PET particles 
are spherical. Accordingly, the surface area of PET N,, changes proportionally 


to Von . The initial PET particle diameter was measured with the help of an 
optical microscope (see Fig. 5.7). In [90], it is assumed that the PET particle 
surface remains constant. Different to [90], a progressive decrease of the 
particles surface is considered in the following model approach. This should 
allow a more precise determination of the reaction kinetics. 


Figure 5.7: Microscope pictures of the PET particles 


The degree of PET conversion a can be calculated from the measured tan 6 
and temperature T (Fig 5.2) with the help of Fig. 5.6. The calculated a can be 
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5 Modelling of the reaction kinetics based on tan 6 


assigned to the corresponding t with the help of Fig 5.3. The time interval 
between two measurements is approximately At - 0.5 s. It is assumed that 
there is no temperature change between two measurements points. 
Subsequently, on the basis of the discretely estimated time- and temperature- 
dependent a, Aa /At can gradually be estimated by the different temperatures. 


A table can be created in which Aa/At, T and N,, data are assigned to each 
other. By solving Eq. 5.7 to k and substituting the variables from the above 
table, the Arrhenius plot can be calculated. The Arrhenius plot shows the 
logarithm of the reaction rate constant Ink over the reciprocal of the 
temperature 1/T [94]. Fig. 5.8 shows the corresponding Arrhenius plot. In this 
figure, the measured data is fitted with a straight line, because, according to 
Eq. 5.5, for the numerical value {x} of k, In{x} is linearly dependent on 1/T. 
The maximum deviation of x as a function of the temperature from the fitted 
results is < 4 %. 
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Figure 5.8: Arrhenius plot for the alkaline hydrolysis of PET 


1 
nfk} = -6014 +5 (5.8) 


The slope of the Arrhenius plot corresponds to - E„/R. From Eq. 5.8 results 
the activation energy E, from PET to 50+5 kJ/mol and the 
Arrhenius constant calculated from the intercept to be 
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5.3 Verification and errors 


exp(5) = 148 + 40 L/(min -cm?). These numbers fit well with the one 
referred in [90]. 


5.3 Verification and errors 


In order to validate the model parameters in section 5.2, three different start 
reaction mixture from Table 4.1 are experimentally examined. 
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Figure 5.9: Measured and modelled dielectric properties for three different start reaction mix- 
tures from Table 4.1 
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5 Modelling of the reaction kinetics based on tan 6 


During the experiment, the heating rate was chosen to be 3 °C/min. At first, 
the numerical solution of Eq. 5.7 is used to model the course of the reaction. 
Subsequently, the tanö course is modelled with the help of Fig. 5.6. The 
corresponding tanö for these three samples are shown in Fig. 5.9. The 


modeled tan 6 progress for the different start reaction mixtures are shown in 
Fig. 5.9. 


From Fig. 5.9 a good agreement between the modelled and measured 
dielectric properties is evident. The most important sources of error are the 
inaccuracy of the measured dielectric properties, the inaccuracy of the modeled 
temperature depended dielectric properties of the unreacted mixture, the 
assumed PET particle thickness and finally the assumption that the losses 
change only due to the variation of the dissolved NaOH concentration in the 
liquid solution. The total error on the predicted tan ô is estimated to be less 
than 10 %. 
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6 Physical and technical design of 
an industrial Wave applicator 


6.1 Concept 
The applicator has to meet the following specifications: 


e It must enable a continuous operation with the use of a metallic 
Archimedean screw within a tube, which must be transparent to micro- 
waves, e.g. glass. 


e It has to provide a homogeneous electric field distribution in axial and 
azimuth direction of the processed material. 


e It has to provide a production rate for the depolymerization of PET of 
at least 12.5 kg/h and accordingly an overall flow rate of at least 
100 kg/h in total for PET, MEG, NaOH and H,0. 


6.2 Archimedean screw geometry 


For the microwave assisted depolymerization process, PET and appropriate 
catalysts have to be transported through the microwave field. This is achieved 
by the use of a metallic Archimedean screw placed in a microwave transparent 
glass tube. To make sure that the metallic structure of the Archimedean screw 
allows an efficient penetration of the electromagnetic field into the liquid, the 
influence of the Archimedean screw design on the microwave power 
absorption has been investigated. For optimization the following geometric 
parameters of the screw were varied: the pitch distance d, and the thread depth 
h,. To simplify the discretization with an orthogonal grid and to reduce the 
simulation efforts the screw has been modelled by a metallic core and parallel 
discs with a distance d,. The simplified screw model is shown in Fig. 6.1. 
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6 Physical and technical design of an industrial u Wave applicator 


Figure 6.1: Archimedean screw (left) and its simplified model for the simulations (right) [95] 


It is assumed that the screw parameters have a comparable influence on the 
total power efficiency, independent of the exact cavity geometry. For this 
investigation a rectangular cavity, similar to that of a household microwave 
oven (310 mm x 375 mm x 200 mm) is used. It is expected that a variation of 
the screw core diameter should have insignificant impact and, therefore, it is 
set equal to 100 mm for the further investigations. In the electromagnetic 
simulation, the screw material is consisted as a perfect electric conductor 
(PEC) and the dielectric constant €, of the surrounding liquid has been chosen 
equal to 30, because this value corresponds to the order of magnitude of 
permittivity of the later processed reaction mixture in cold state. The 
corresponding wavelength A at the magnetron frequency of 2.46 GHz inside 
the material is equal to 22.3 mm. The depth of the thread h, as well as the pitch 
d, were varied between 1/4 and A in order to study both the resonant and non- 
resonant behavior. Note that for a fixed external screw diameter, the material 
volume depends on the final geometry of the screw. 


The estimated absorbed power in the liquid material (given in percentage of 
the input microwave power) for €’, = 30 and tan ô = 0.01 or tan ô = 1 are 
shown in Table 6.1 for different d, and h,. Both parameters are scaled with the 
wavelength A = 22.3 mm. The absorbed power is integrated over the volume 
of the material. 


Table 6.1 reveals a significant influence of the screw geometry on the 
efficiency of the applicator. In particular for a material with a small dielectric 
loss constant (tan ô = 0.01) the screw depth h, of A/4 and 3 A/4 ‚as well as 
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6.2 Archimedean screw geometry 


an increasing pitch d,, enables a higher absorbed power. Furthermore, 
Table 6.1 shows that the geometry parameters of the screw are much more 
important at lower loss tangents, due to the higher penetration depth of the 
microwaves. The dielectric loss tangent of the liquid used in the process is in 
the range of 0.3 <tan ô < 1, and correspondingly the differences in the absorbed 
power between resonant (grey lines in Table 6.1) and non-resonant screw 
parameters are no longer obvious and may almost disappear for tan ô = 1, if 
the penetration depth D, is in the range of about 3 mm. 


Table 6.1: Simulated relative power absorption for varius screw parameters, scaled to A for a 
dielectric material with €, = 30 and tan ô = 0.01 and 1 respectively [95] 


Dimensions Pabs [%] 

d, hs tan ô = 0.01 tand = 1 
AA AA 11 84 
a4 2 6 84 
a4 31/4 24 84 
M4 À 15 84 
A/2 AA 40 75 
M2 /2 20 76 
r/2 30/4 41 77 
2 À 40 77 
31/4 1/4 41 84 
31/4 2 34 82 
31/4 30/4 66 82 
31/4 x 72 82 

À a4 92 86 

À A/2 24 82 

À 31/4 86 83 

À À 45 83 
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6 Physical and technical design of an industrial u Wave applicator 


It is found that a screw pitch of d, = n: 4/2; n = 1,2,3 in combination with 
a depth of the thread of 


A 
h, = (2n +1): qi n= 012,.. (6.1) 


enable the best efficiency. This can be understood from the fact that the 
tangential component of electric field at a perfect conductive surface is equal 
to zero and has a maximum at a distance of h,. The screw acts like a 
quarter-wave wavelength transmission line, equivalent to a short to open 
transformation and enables a maximum coupling of the microwave power into 
the material in-between two neighboring threads. 


Figure 6.2 shows the relative absorbed power for different tand for an 
optimized and a completely unfavorable screw geometry. It is obvious that the 
improper selection of the geometrical parameters of the Archimedes screw can 
lead to a significant reduction of the overall system efficiency. This factor is 
especially important for tan ô < 0.5. This optimization is proposed in [95] 
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Figure 6.2: Simulated relative power absorption as a function of 
tan ô at 2.46 GHz and €, = 30 [95] 
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6.3 Energy required for a specific volume 


6.3 Energy required for a specific volume 


The reaction temperature ofthe depolymerization of PET lies in the range from 
150 °C to 220 °C [28]. The specific heat capacity of the reaction mixture Cp ges 


can be estimated as weighted average over the individual reactants as follows 


= MpeTCpper + MmEsCpyes + ™naoH"pnaon + M0 pH,0 (6.2) 


C. = 
Pges 
g Mpgr + Muec + Myaon + Mio 


The specific heat capacities and the quantities of the individual reactants are 
shown in Table 6.2. 


Table 6.2: Specific heat capacity of the various reactants and the corresponding mass in the 
used reaction mixture. 


Material Cy [kJ/(kg K)] mass [g] 
PET 1.1 [96] 180 
MEG 2.2 [74] 955 
H20 4.2 [74] 140 

NaOH 1.49 [97] 82.4 


Accordingly, the mentioned reaction mixture provides a specific heat capacity 
of approximately 2.2 kJ/(kg K). The power Pheat required to raise the 
temperature by AT of material with a total mass m within the time At can be 
calculated as follow [98]: 


Preat = Cpges MAT /At (6.3) 


In the application considered here, with a reaction mixture mass of 100 kg and 
a start temperature of 20 °C, the required power is between 8.0 kW and 
12.3 kW depending on the final reaction temperature, without consideration of 
thermal heat losses. 
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6.4 Estimation of the applicator segments 


The outer diameter of the glass tube is chosen to be 250 mm with a wall 
thickness of 10 mm. The mentioned glass diameter has been selected according 
to the maximum dimensions of commercial quartz glass tubes approved for 
pressure application. The modular concept of the reactor implies a subdivision 
of the applicator’s geometry into N equivalent segments. With the use of this 
kind of modular concept and the limitations of the glass tube, any specific 
temperature profile along the reactor may be realized. It requires the control of 
microwave power level in each segment individually, thus resulting in the 
corresponding temperature scenario along the material axis (see Fig. 6.3). 


Segment 


Z 


Figure 6.3: Modular design of the microwave reactor (top), temperature profile along 
the MUT (bottom) [95] 
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In Fig. 6.3 the temperatures T,, Tz, ... Ty are the target temperatures at the end 
of each segment. In the particular reaction the reaction mixture should be 
heated up to the reaction temperature as fast as possible and then the 
temperature should be kept constant until the reaction is complete. 
Accordingly, two segments would be sufficient. 


Another crucial criterion for the number of segments required, is the size of 
the available microwave generators (magnetrons) and their prices. As the 
commercial available magnetrons with more than 3 kW nominal microwave 
power are equipped with an additional circulator, which is very expensive, only 
magnetrons with a nominal power from | kW, 2 kW and 3 kW are considered. 
In section 6.3 the power required to bring the reaction mixture to the reaction 
temperature is estimated between 8 kW and 12 kW. As the magnetrons with a 
nominal power between 2 kW and 3 kW are the most economical in terms of 
price per watt. Four segments would be preferred. 


Since the process temperature exceeds the boiling point of the aqueous phase 
of the reactive mixture, a pressurized microwave reactor is required. This 
reduces the failure probability of the glass tube. Thus, the preferable geometry 
of the microwave cavity is cylindrical. At the same time, considering the 
investment and installation efforts, the number of coupling ports with pressure 
sealed microwave windows and accordingly the number of microwave 
generators must be reduced to the minimum. Therefore, in the present design 
only one magnetron and hence one quartz windows per segment is considered. 


6.5 Applicator design 


In this section, the design of the industrial scale microwave cavity based on the 
measured material properties is presented. Fig. 6.4 shows the simulated electric 
field density in a cylindrical cavity with 336 mm in diameter and a length of 
250 mm with a single microwave port on the left side. The core of the 
Archimedean screw occupies the hollow space in the cavity center. The 
cylindrical sample is enclosed by the glass tube with a diameter of 250 mm. 
From Fig. 6.4 it can be seen, that the most of the microwave field is absorbed 
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through the material near the port position and that the electric field shows a 
strong inhomogeneity in the z-direction as well as in the azimuthal direction. 


Standard waveguide (WR430) Sample 


| 


Beeessses 


< 


250 mm x» 


Figure 6.4: Simulated electric field in the cylindrical cavity [95] 


The inhomogeneity in the z-direction can be avoided by limiting the cavity 
length in z-direction to 50 mm < A/2. In this way a constant E-field is forced 
in z- direction. The dominant mode in that case is a TE,9g, -mode. The 
corresponding simulation results are shown in Fig, 6.5. 
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[v/m (log) 1 


Z 
IE B Le x 
Figure 6.5: Electric field density in a ‘short’cylindrical microwave cavity with electric length 
in z-direction < A/2 [95] 


It should be noted, that in the ‘short’ cavity the azimuthal inhomogeneity of 
the electric field remains. A further improvement of the field distribution is 
realized through an optimized power distribution network design. Such an 
approach enables a well-defined power distribution over the azimuthal 
direction of the resonator. The optimized distribution network is shown in the 
Fig.6.6. 
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Matching 
element 


coupling 
ports 


Figure 6.6: Simulated geometry of the power distrubution network: center zone (left), 
side view (right) [95] 


This design consists of a single microwave port, which is coupled to the 
external resonator circuit. The outer diameter of the resonator is 534 mm. The 
resonator structure has a total effective length of 250 mm and subdivided into 
5 zones of equal length (Fig 6.6). Furthermore, looking in the azimuthal 
direction of each zone four equidistant coupling ports are designed to couple 
the microwave power to the corresponding heating cavity. It is particularly 
important that the phase velocity between the actual heating cavity and the 
external resonator is the same. The number of zones of equal size in a single 
segment can be varied depending on the required sample length and/or required 
power density. In that case the coupling apertures between the inner and the 
outer resonator have to be modified in an appropriate way. The coupling 
apertures between the different zones should have a width of about A and are 
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only depending on the width of the outer resonator. The width and length of 
this coupling aperture for an outer diameter of the resonator of 534 mm are 72 
mm and 157 mm, respectively. In the middle zone this coupling is 
positioned opposite ofthe waveguide port (see Fig. 6.6 left). And for the other 
zones, it is positioned axially symmetrical. 


Figure 6.7: Electric field distribution of the power distribution network design in azimuthal 
direction of zone c (upper right); of levels b and d (lower left); of levels a and e 
(lower right) and in axial direction (upper left) [95] 


Table 6.3 shows the power distribution at the microwave cavity for the five 
individual segments (a, b, c, d, e). The coordinates of the different radial MUT 
sections is shown in Fig. 6.7. The listed power is given in percent of the total 
power absorbed in the applicator. 
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Table 6.3: Power distribution in the MUT located in the power distribution network 
[% ofthe absorbed power] [95] 


MUT section a b c d e 

(coordinates) 
1 2.5 2.4 2.5 2.4 2.5 
2 2.4 2.6 2.5 2.6 2.4 
3 2.2 2.8 2.4 2.8 22: 
4 2.3 2.9 2.5 2.9 2.3 
5 2.2 2.8 2.4 2.8 2.2 
6 2.4 2.6 2:5 2.6 2.4 
7 2.5 2.4 2.5 2.4 2.5 
8 2.8 2.2 27 2.2 2.8 


The proposed design shows some similarities with the one presented from [22]. 
In [22] a so-called ‘Coaxial Wrap-Around TE,, Mode Converter’ is 
introduced. This is based on two coaxially arranged WR90 waveguides that 
shares a common broad-wall. The design is very similar to the middle segment 
(segment c) of Fig. 6.7. However, there are a number of key differences with 
the concept presented here. One of them is the organizing of the coupling 
between the inner and outer waveguide. In the here proposed design, it is 
realized via a sidewall of the waveguide contrary to [22] where the coupling is 
organized through a broadwall. This enables a fully constant electric field in 
the axial direction. Moreover, the coupling slots between the different cavity 
segments are tapered to enable a well-defined power distribution. Last but not 
least the proposed designed consists of many ‘stacked’ zones. This enables a 
constant and defined field distribution over any material length. It is not limited 
to samples with a length < 2/2. 


The proposed design enables a good and well-defined field distribution (see 
Table 6.3). Since the transition from the waveguide port to the ring-shaped 
external resonator acts like a waveguide power divider, an additional matching 
element is needed to achieve a maximum system efficiency. The implemented 
optimized matching element in the front of the waveguide port (only in the 
middle zone) is shown in Fig. 6.6. It must be oriented vertically to the port. 
Already a small misalignment (<3°) of this element leads to a significant 
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change of the power distribution in the middle zone c. The power distribution 
in adjacent zones is not affected. For example, a rotation clockwise of the 
matching element from 1° to 3° leads to a difference in the absorbed power 
between sample section 1 & 7 by factor of 1.5 to 3 respectively. The total 
absorbed power in the applicator is reduced by 2 %. This matching element 
leads to significant reduction of the reflected power. The reflected power with 
and without “matching element’ is presented in Fig. 6.8. 


— BEE - — 


without matching 
[with matching 


root} 


-20 


-25 i L i L fi 
2.4 2.41 2.42 2.43 2.44 2.45 2.46 2.47 2.48 2.49 2.5 
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Figure 6.8: Simulated reflected power at the input of the power distribution network with and 
without matching element versus frequency 


Fig. 6.9 shows the total absorbed power versus frequency for different 
dielectric constants of the material. The design is characterized by a high 
robustness against material permittivity changes (verified for high loss 
materials with tan ô > 0.3). This can be explained by the fact that the 
resonance frequency of the complete resonator is given by the external power 
distribution network. In the actual heating cavity is existed a running wave. 
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Pabs [%] 


2.4 2.45 2.5 


Frequency [GHz] 


Figure 6.9: Absorbed power for different e/ and tan ô = 0.5 versus frequency 


To avoid a physical electric contact between the cylindrical metal plates 
(levels separators) and the outer metal tube, a 4/4 = 47 mm stub filter may be 


used. This filter makes a virtual open -to -short transformation. Here the gap 
between filter and metal tube is < 2 mm. 


>l le 


Figure 6.10: 4/4 filter structure to enable a contactless operatic _ 
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Pabs [%] 


2.4 2.45 2:5 


Frequency [GHz] 


Figure 6.11: Absorbed power versus frequency for different gap thicknesses of the 
A/4 filter structure 


At larger distance between inner structure and outer metal tube, the operating 
mode is strongly suppressed and the matching between the source and the load 
is lost. This results in a significant loss of system efficiency as shown in the 
following figure for a 2 and 4 mm gap, respectively. As an alternative contact 
springs may be used although those are not specified for such high-power 
applications [99]. Depending on the spring geometry contact gaps in the range 
from 0.14 mm to 0.7 mm could be bridged. 


6.6 Performance of the fabricated applicator 


6.6.1 Test set 


The design presented in section 6.5 is based on a complex geometry, which 
consist of several thin sheets of metal that need good electrical contact with 
each other. Because of the complexity of the geometry, a prototype is very 
costly to manufacture. A prototype production based on 3d metal printing is 
out of the question due to the excessive size of the cavity. To enable a 
cost-efficient manufacturing the design of the prototype has been slightly 
modified. Modifications at the following elements are done: the wall 
thicknesses and radius of the edges. For testing purposes, a carbon blended 
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silicone rubber from Gummiwerk Kraiburg Gmbh & Co. KG [100] is chosen. 
The thermal conductivity of the chosen MUT is low to prove the field 
homogeneity. A hollow cylindrical sample with 183 mm outer diameter and 
171 mm inner diameter and dielectric parameters of €’ = 15,tanö = 0.3 is 
chosen. The MUT consists of several stacked cylinders (layers) with tiny 
vertical cuts. Hence, it forms a sandwich structure. With the vacuum bag the 
MUT is fixed on the metal core (sample holder, see photo on the Fig. 6.12). To 
check the homogeneity of the absorbed microwave power in the load, the 
surface temperature of the sample is recorded with an infrared camera. Firstly, 
the sample is exposed to the microwave inside the designed applicator. In the 
following step, it is removed from the applicator and a thermal imaging is 
recorded. For the experimental verification a magnetron with an output power 
of 800 Watt and the infrared camera FLIR type AX5 are used. Sample pictures 
of the manufactured microwave cavity are shown in Fig. 6.12. 


Figure 6.12: Fabricated microwave cavity [95]. 
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6.6.2 Results 


For the verification of the developed applicator, the final CAD file from the 
manufactured applicator is simulated with CST MWS. The comparison with 
the measurement is shown in Fig. 6.13. From Fig. 6.13 is obvious that the 
resonance of the developed applicator match very good to the simulation 
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Figure 6.13: Simulated and measured S-parameters 


Fig. 6.14 compares the measured temperature distribution, against the 
simulated absorbed microwave power distribution within the sample. The 
measured temperature inhomogeneity is found to be +8°C at a mean 
temperature level of 90 °C (see Fig. 6.14). A proper agreement between 
simulation and experiment exists. It must be noted, that in the experiment the 
sample does not completely fill the glass tube inside the applicator. It enables 
a fast temperature measurement but results in some differences with the 
simulated absorbed power profile. 
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Figure 6.14: Simulated power distribution (left), measured temperature distribution (right) [95] 


6.7 Process optimization 


The goal is to find the optimum power level of the power sources individually 
for each segment to get a full-reacted mixture at the end of the applicator. 
Based on the simulated power distribution in a single segment (see Table 6.3), 
the demanding temperature profile of the sample passing the microwave cavity 
is estimated for a four-segment configuration with a total length of one meter. 
The capacity of the proposed reactor is fixed to about 100 kg/h. Accordingly, 
the Archimedean screw angular and axial speed is 0.067 rotations per second 
and 2 mm per second. For the optimization, it is assumed that the thermal 
convective losses of the mixture remain homogeneous in each zone. This 
assumption is feasible because the applicator consist ofa total of 20 zones (four 
segments with five zones each). The MUT is heated over its entire length by a 
AT < 180 °C. Hence, it is expected a temperature deviation of about +4.5 °C 
within a zone. 


As magnetrons are considered as power sources. The following levels of 
nominal power are industrially available: 1, 2, 3, 6 and 10 kW. Ideally, the 
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maximum power of the respective power source is used. The optimization 
bases on the reactions kinetics developed in section 5 and results ina demanded 
power of 6 kW in the first segment, 2 KW in the second segment and 1 kW in 
the third and fourth segment. The convection losses through the glass tube over 
the entire applicator can be determined using Eq. 3.9. And for this temperature 
scenario, they amount to approx. 660 W in total. 
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Figure 6.15: Temperature profile and reaction progress across the applicator at a reaction tem- 
perature of < 180 °C 


In the above scenario, the reaction temperature is limited to 180 °C. Under this 
assumption, the minimum reaction time is around 5 minutes and 30 seconds, 
excluded any heating up time. Therefore, a further increase of the installed 
power does not follow to significantly higher yields. Since the reaction may 
take place at temperatures up to 220 °C, it is theoretical feasible to double the 
screw speed and accordingly the reactor capacity by increasing the installed 
power levels of 10 kW in the first and second segment, 3 kW in the third 
segment and 0.3 kW in the fourth segment. This scenario follows to very 
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critical reaction pressure conditions. The corresponding reaction profile is 
shown in the following figure. 
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Figure 6.16: Temperature profile and reaction progress across the applicator at 
temperatures < 220 °C 
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Polyethylene terephthalate, PET, is a widely used plastic. It shows excellent 
thermal and mechanical properties and chemical resistance. On the one hand, 
this robustness of PET is desirable, but on the other hand it is also detrimental 
for the environment if it is disposed of after use. This makes PET recycling 
indispensable. PET recycling can be achieved by mechanical processing, 
chemical recycling and energy recovery. Through chemical recycling, the 
molecules that originally make up the polymer (PET) can be completely 
recovered without degradation. This technique is not widely used due to the 
long reaction time and the high-energy demand. The microwave-assisted 
depolymerization of PET is expected to take place in significantly shorter 
reaction times and under milder reaction conditions (e.g. pressure 
and temperature) [12]. This is probably because in case of microwave-assisted 
depolymerization, apart from a temperature effect, the microwaves also initiate 
currents in the material, which leads to a higher mobility ofthe Na* ions. 


In this thesis, an innovative microwave applicator for PET recycling is 
developed with the aim to build an efficient, large-scale process that allows 
continuous operation. The development is carried out for the Industrial, 
Scientific and Medical (ISM) frequency band at 2.45 GHz. 


The starting point for the successful development of the microwave applicator 
is the knowledge of the dielectric properties of the materials involved at the 
process frequency and temperatures. Therefore, an advanced, fully automated 
system for accurate in-situ monitoring of dielectric properties during 
microwave heating has been developed. A dual-mode approach is combined 
with a 3D full-wave simulation. The system allows measurements in wide 
ranges of the dielectric constant 15 <e,<80 and loss tangent 
0.01 <tanö < 1. The coupling design allows an isolation of about 60 dB 
between both operating modes, thus the interference can be neglected. 
Furthermore, the sample holder is sealed at both ends with connectors that 
withstand a pressure of 18 bar. This allows the dielectric characterization of 
liquids above the boiling point. 
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The measurement of the time- and temperature-dependent dielectric loss 
tangent enables the determination of the actual reaction course. Based on these 
data a model for the reaction kinetics is developed. The model allows an 
optimization of the process parameters, such as heating time and temperature 
holding time. This model also allows the definition of the power requirement 
per reaction segment to enable a targeted material throughout under 
consideration of the convection losses across the glass tube. 


Accurate power measurements enable a qualitative calorimetric estimation of 
the PET depolymerization process. Based on RF and multi-physics modeling 
of the experimental setup, the specific heat capacity of well-known liquids is 
determined based on a microwave calorimeter for the first time. 


An advanced microwave cavity on an industrial scale, which enables power 
efficient microwave heating of liquids with high permittivity at 2.45 GHz, has 
been developed. The developed applicator is modular and allows 
homogeneous heating of a large cylindrical sample with an outer diameter of 
250 mm. The sample volume used in the simulation is 11 liter. The use of a 
coupled resonator structure allows a specific power distribution in the cavity. 
In addition, a matching element at the port position leads to a significantly 
increased power efficiency. 


The reached power efficiency is > 95 % for high permittivity materials in the 
range from 15 < e€, <25 and tanö > 0.3. The maximum temperature 
deviation is estimated to be below 4 °C for a sample volume of 11 liter. 


Due to its modular design the applicator can easily be adapted to the 
requirements of different reaction conditions and different reactants. The first 
experiment shows a very promising homogeneity of the absorbed RF power 
(+8 °C). The temperature distribution corresponds well to the expectations 
from the simulation for the MUT geometry used. The continuous operation of 
the applicator is realized with an Archimedean screw, which transports the 


displacement through the microwave pattern. The screw geometry was 
investigated and optimized independently of the final design of the applicator. 
It is shown that the screw geometry becomes especially important for materials 
with tand < 0.3. 
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Based on the work done in this thesis, the calorimetric measurements setup 
could be further optimized in order to improve the power measurements 
accuracy. For example, a more stable bidirectional coupler and or improved 
manufactured cavity, which allows a more precise definition ofthe wall losses, 
would significantly improve the accuracy of the calculated absorbed power. 
Furthermore, a faster temperature sensor would allow faster heating ramps and 
more accurate temperature control. 


Regarding the industrial scale applicator, even lower frequencies, such as 
915 MHz, would be interesting. This enables higher penetration depth of the 
microwaves into the reaction mixture and subsequently allows the core of the 
Archimedean screw to be reduced. In this way, a significant increase of the 
applicator capacity tube may be reached by use ofthe same glass tube diameter. 


In the frame ofthe European Union’s Horizon 2020-SPIRE-2015 research and 
innovation program SYMBIOPTIMA in cooperation with the industry, a 
microwave applicator similar to the one presented in this thesis was also 
developed and is successfully installed at one of the industry partners. They 
have shown its efficiency and the simulations match very well. Based on the 
results of this applicator, one of the project partners won the 2018 Innovation 
Radar Prize [101] 
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Frequency and temperature dependent dielectric properties of MEG 
and NaOH 


Fig 8.1 shows the measured dielectric parameters of MEG for three different 
temperatures. It can be seen that the dielectric constant €, of MEG increases 
with the temperature and the dielectric losses decreases (see Fig 8.1). 
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Figure 8.1: Measured dielectric properties of MEG using the coaxial probe HP 85070B at 
three different temperatures in a wide frequency range 
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The Debye model can describe the behavior ofthe measured permittivity from 
Fig 8.1. Using the Debye model, the frequency dependent permittivity can be 
determined as follows [102] 


4 fs fm 
os (1 + jot)? (1) 


i: 


E€ = 


_ WT(Es — Eco) 
~ (14+ jwt)? 


A 


(8.2) 


where €,.corresponds to the permittivity of the medium at very high 
frequencies, €, the low-frequency permittivity of the medium and r the 
characteristic relaxation time of the medium. T is connected with the 
temperature as follows: 


Eq 
T= To ExPP) (8.3) 


where To is the relaxation constant, Eg is the activation energy and kp is the 
Boltzmann constant. For T = w"! the course of e’follows a step from €, to E€% 
and €” runs through a maximum at the same frequency. From eq. 8.3, it is clear 
that an increase in T results in a shift of the frequency dependent permittivity 
to higher frequencies (to the right). This behavior follows the curves measured 
here (see Fig 8.1). 


Fig 8.2 shows the measured dielectric properties of a 5 wt.% NaOH solution 
in water. It can be seen that the dielectric losses increases with the temperature. 
It should be noted that the losses in the mixture are mainly dominated by the 
conductivity losses, which normally increases exponentially with temperature. 
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Figure 8.2: Measured dielectric properties of 5 wt.% NaOH with H2O mixture using 
the coaxial probe HP 85070B at three different temperatures in a wide 
frequency range 


Modelling of the reaction kinetic 


A usage of the curve-fitting tool from MatLab ‘cftool’ shows that the behavior 
of the unreacted reaction mixture over the temperature can be approximated by 
a Gaussian function with five terms. The expected error on the loss factor due 
the fitting is estimated to about 0.11 %. 
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5 


(2x)? 
tan Östart (T) = axe “x (8.4) 


x=1 


Table 8.1: Course of the coefficients and their approximation from the unreacted 
reaction mixture 

x ax bx Cx 

1 0.1980 15.1481 2.7032 

2 0.3301 14.1132 16.1049 

3 0.1642 40.8480 30.7718 

4 0.0227 75.7096 13.8588 

5 0.6254 141.1441 430.6950 

5 
Se 
tan Ögna (T) = aye \ & (8.5) 
x=1 


The course of the coefficients and their approximation for the reaction product 
are given in Table 8.2 and from the reaction mixture during the reaction in 
Table 8.3. The expected errors through this fitting for the reaction product and 
the reaction mixture are estimated to be about 0.18 % and 0.1 % 
correspondingly. 
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Table 8.2: Course of the coefficients and their approximation from the reaction product. 


x ax bx Cx 

1 0.0475 27.6109 7.4211 
2 0.3781 207.5609 106.1568 
3 0.0109 123.8985 13.4328 
4 6129.5 -325.851 109.2479 
5 0.6894 -637.99 633.3769 


Table 8.3: Course of the coefficients and their approximation from the loss factor during the 


reaction. 
X ax bx Cx 
1 41.4544 -30.5646 21.5679 
2 0.1496 26.3432 6.8576 
3 0.6486 -4.5765 67.3126 
4 0.0153 183.9751 14.9851 
5 0.4651 178.0238 152.0581 


The fitting of the time depending loss factor during the reaction is realized also 
with the help of a Gaussian function with 8 terms. The expected error through 
the fitting is 0.15 % 


8 


; | t-by\? 
tandreaction(t) = eax Cx ) (8.6) 


x=1 
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Table 8.4: Course of the coefficients and their approximation from the time depending loss 
factor of the reaction mixture. 

x ax bx cx 

1 3.798ell -4289.3 825.3049 
2 0.3269 196.8953 159.6712 
3 0.2332 400.6403 129.1403 
4 0.0906 621.5689 256.4009 
5 0.4440 3492.2 1046.7 
6 0.4824 752.6186 1623.7 
7 0.0236 2710.1 419.3913 
8 0.1099 2258.3 606.1110 
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By using/microwave heating, the alkaline hydrolysis 
polyethýlene terephthalate (PET) can be significantly ac- 

celerated from about 3 hours to less than minutes, 
while the reactiön.takes place at lower fepipsratures un- 
der microwave conditions. However, the benefits have 
only/ been demonstrated’ a laboratory scale to date. 
This has led to the need to develop a suitable industrial 
migrowave reactor to enable indu trial application of this 
reaction. For the first time, this work. resents an novel 
industrial microwave appli ator at 2. ash with homo- 
geneous distribution to support the alkaline hydrolysis of 
PET, which allows an efficient and continuous operation. 
In addition, an innovative dielectric and calorimetric meas- 
urements setup is presented. Furthermore, the modelling 
of the reaction kinetics based on the measured die hectric 
parameters is presa is N 
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